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We investigate properties of an energetic atom propagating through strongly interacting atomic 
gases. The operator product expansion is used to systematically compute a quasiparticle energy and 
its scattering rate both in a spin- 1/2 Fermi gas and in a spinless Bose gas. Reasonable agreement 
with recent quantum Monte Carlo simulations even at a relatively small momentum k/kp > 1.5 
indicates that our large-momentum expansions are valid in a wide range of momentum. We also 
study a differential scattering rate when a probe atom is shot into atomic gases. Because the 
number density and current density of the target atomic gas contribute to the forward scattering 
only, its contact density (measure of short-range pair correlation) gives the leading contribution to 
the backward scattering. Therefore, such an experiment can be used to measure the contact density 
and thus provides a new local probe of strongly interacting atomic gases. 

PACS numbers: 03.75.Ss, 03.75.Nt, 34.50.-s, 31.15.-p 



CONTENTS 

I. Introduction 1 

II. Summary of results and discussions 

A. Quasiparticle energy and scattering rate 

B. Differential scattering rate 

C. Comparison with Monte Carlo simulations J 

III. Spin- 1/2 Fermi gas 

A. Formulation H 

B. Operator product expansion 3 

C. Wilson coefficients 8| 

D. Expectation values of local operators H 

E. Single-particle Green's function ^ 

F. Three-body problem [It] 

IV. Spinless Bose gas 13 

A. Formulation 12| 

B. Operator product expansion 12| 

C. Wilson coefficients 13| 

D. Single-particle Green's function 13| 

E. Three-body problem 

V. Differential scattering rate HI 

A. Spin- 1/2 Fermi gas HI 

B. Three-body problem 16| 

C. Differential scattering rate Ij] 

1. Contribution of number density 13 

2. Contribution of current density 18| 

3. Contribution of contact density HI 

D. Spinless Bose gas 2t] 

E. Three-body problem [2l| 



VI. Weak-probe limit 

A. Spin- 1/2 Fermi gas 

B. Dynamic structure factor 

C. Differential scattering rate 

D. Spinless Bose gas 



VII. Conclusions 

Acknowledgments 

A. Derivation of Wilson coefficients 

1. One-body and two-body sectors 

2. Three-body sector 

B. Details of solving integral equations 

1. Spin- 1/2 Fermi gas 

2. Spinless Bose gas 



25 
2G 



28 
28 
30 



[31| 



C. Derivation of optical theorem in Eq. (|5.24[) 
References 

I. INTRODUCTION 



Strongly interacting many-body systems appear in var- 
ious subfields of physics ranging from atomic physics to 
condensed matter physics to nuclear and particle physics. 
An understanding of their properties is always impor- 
tant and challenging. Among others, ultracold atoms of- 
fer ideal grounds to develop our understanding of many- 
body physics because we can control their interaction 
strength, dimensionality of space, and quantum statis- 
tics at will [U-Il]. Many-body properties of strongly in- 
teracting atomic gases have been probed by a number 
of experimental methods including hydrodynamic expan- 
sions j Bragg spectroscopies |9Hl2T] , radio- frequency 
spectroscopies [13H15j . precise thermodynamic measure- 
ments [l6l422| . and collisions of two atomic clouds (23T - 
[25j . In particular, the photoemission spectroscopy em- 
ployed in Refs. [2g, [27} is a direct analog of that known 
to be powerful in condensed matter physics [28| . 

On the other hand, often in nuclear and particle 
physics, high-energy particles play important roles to re- 
veal the nature of target systems. For example, neutron- 
deuteron or proton-deuteron scatterings at intermedi- 
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FIG. 1. Schematic of a proposed scattering experiment in which we shoot a probe atom \ into an atomic gas with a large 
momentum k and measure its differential scattering rate. Three leading contributions come from two-body scatterings propor- 
tional to the number density n and current density j of the target atomic gas and from a three-body scattering proportional to 
its contact density C (measure of short-range pair correlation). We will find that the number and current densities contribute 
to the forward scattering [6 < 90°) only and therefore the contact density gives the dominant contribution to the backward 
scattering (6 > 90°). 



ate or higher energies are important to reveal the ex- 
istence of three-nucleon forces in nuclei [29| - |3]| . Also, 
two-nucleon knockout reactions by high-energy protons 
or electrons have been employed to reveal short-range 
pair correlations in nuclei [32], |33| . Furthermore, the 
discovery of "jet quenching" (i.e., significant energy loss 
of high-energy quarks and gluons) at Relativistic Heavy 
Ion ColliderJRHIC) [M H| and Large Hadron Collider 
(LHC) [H, H3] is one of the most striking pieces of evi- 
dence that the matter created there is a strongly interact- 
ing quark-gluon plasma [38j |. Such high-energy degrees 
of freedom to probe the nature of a quark-gluon plasma 
are generally referred to as "hard probes" [39j. In con- 
densed matter physics, the use of high-energy neutrons 
to probe the momentum distribution of helium atoms in 
liquid helium was first suggested in Re fs. I4 fj| .l4ll| and has 
been widely employed in experiments [42j-|44|. 

Now in ultracold-atom experiments, analogs of the 
hard probe naturally exist because a recombination of 
three atoms into a two-body bound state (dimer) pro- 
duces an atom-dimer "dijet" propagating through the 
medium. Such three-body recombinations rarely occur 
in spin- 1/2 Fermi gases because of the Pauli exclusion 
principle [45[ but frequently occur in spinless Bose gases 
and have been used as a probe of the Efimov effect [46l — 
[50j . While atom-dimer "dijets" are normally consid- 
ered to simply escape from the system, multiple colli- 
sions of the produced energetic dimer with atoms in the 
atomic gas were argued in Ref. [ItJ to account for the 
observed enhancement in atom loss. A well-founded un- 
derstanding of collisional properties of an energetic atom 
or dimer in the medium is clearly desired here (see early 
works [5l], H2] and also recent one [EH). In addition to 
these naturally produced energetic atoms, it is also possi- 
ble to externally shoot energetic atoms into an atomic gas 
in a controlled way and measure the momentum distri- 
bution of scattered atoms. Indeed, closely related exper- 
iments to collide two atomic clouds have been performed 
successfully for Bose gases [54146 lj and strongly interact- 
ing Fermi gases [23l425l |. which have been analyzed theo- 
retically [62H65| . 

In this paper, we investigate various properties of an 



energetic atom propagating through strongly interacting 
atomic gases. Such properties include a quasiparticle en- 
ergy and a rate at which the atom is scattered in the 
medium. Both the quantities reflect many-body proper- 
ties of the atomic gas and, in particular, the scattering 
rate may be useful to better understand multiple-atom 
loss mechanisms due to atom-dimer "dijets" produced by 
three-body recombination events (47l . [53| . Also we pro- 
pose a scattering experiment in which we shoot a probe 
atom into the atomic gas with a large momentum and 
measure its differential scattering rate (see Fig. [TJ. Re- 
sulting scattering data must bring out some information 
about the target atomic gas. What can we learn about 
the strongly interacting atomic gas from these scattering 
data? This question will be addressed in this paper. 

Seemingly, these problems are difficult to tackle be- 
cause of the nature of strong interactions. Quite remark- 
ably, however, these problems can be addressed in a sys- 
tematic way. This is because the atom with a large mo- 
mentum probes a short distance at which it finds only a 
few atoms. Therefore, apart from probabilities of finding 
such few atoms in the medium, our problem reduces to 
few-body scattering problems. At the end, we will find 
that such an energetic atom can be useful to locally probe 
many-body aspects of strongly interacting atomic gases. 

Atomic gases created in laboratories at ultralow tem- 
peratures and quark-gluon plasmas created at RHIC and 
LHC at ultrahigh temperatures are both strongly inter- 
acting many-body systems. In spite of the fact that they 
are at two extremes, various analogies have been dis- 
cussed in the literature such as hydrodynamic behaviors 
and small shear viscosity to entropy density ratios [o6| . 
This work intends to build a new bridge between them 
from the perspective of "hard probes." 

Since this paper turns out to be long, we first sum- 
marize our main results in Sec. [Hi discuss their conse- 
quences, and compare them with recent quantum Monte 
Carlo simulations. Our main results consist of the quasi- 
particle energy and scattering rate of an energetic atom 
in a spin-1/2 Fermi gas (Sec. IIII[) . those in a spinless 
Bose gas (Sec. IIV[) . and the differential scattering rate 
of a different spin state of atoms shot into a spin-1/2 
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Fermi gas or a spinless Bose gas (Sec. [V}. Furthermore, 
a connection of our hard-probe formula derived in Sec. [V] 
with dynamic structure factors in the weak-probe limit 
is elucidated in Sec. EU Finally, Sec. I VIII is devoted to 
conclusions of this paper and some details of calculations 
are presented in Appendices lA"rfCl 

Throughout this paper, we set K = 1, /cb = 1, and 
use shorthand notations (k) = (ho, fc), (x) = (t,x), kx = 

kot kx, and tj}' dip = [i/j^ (dip) — (dip <! )ip)/2. Also, note 
that implicit sums over repeated spin indices a = J, J, are 
not assumed in this paper. 



II. SUMMARY OF RESULTS AND 
DISCUSSIONS 

We suppose that atoms interact with each other by a 
short-range potential and its potential range ro is much 
smaller than other length scales in the atomic gas such 
as an s-wave scattering length a, a mean interparticle 
distance n -1 / 3 , and a thermal de Broglie wavelength 
At ~ 1/V mT. Furthermore, we suppose that a wave- 
length of the energetic atom |fc| _1 is much smaller than 
the latter length scales but still much larger than the 
potential range. Therefore, the following hierarchy is as- 
sumed in the length scales: 

(2.1) 



r < Ifej" 1 < \a\, n- 1 / 3 , X 7 



Since the potential range is much smaller than all other 
length scales, we can take the zero-range limit ro — » 0. 
Then physical observables of our interest are expanded in 
terms of small quantities l/(a|fe|), n l ' 3 /\k\, 1 / (Ar|fe|) <C 
1, which are collectively denoted by 0(k~ 1 ), and various 
contributions are organized systematically according to 
their inverse powers of \k\. 



A. Quasiparticle energy and scattering rate 

In the case of a spin- 1/2 Fermi gas with equal masses 
m = m-j- = TO^, the quasiparticle energy and scattering 
rate of spin-up fermions have the following systematic 
expansions in the large-momentum limit (see Sec. Mil for 
details) 



E t (k) 
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1 The energy is often measured with respect to a chemical poten- 
tial. In this case, the quasiparticle energies in Eq. and 1 12.41 
should be replaced with Ef(k) — fif and Ei,(k) — fn, respectively. 



Here, a/ is an s-wave scattering length between spin- up 
and -down fermions, nj, is a number density of spin-down 
fermions, and Cf is a contact density which measures the 
probability of finding spin-up and -down fermions close 
to each other |67H69l | . The results for spin-down fermions 
are obtained simply by exchanging spin indices f ■<-> J,. 

On the other hand, in the case of a spinless Bose gas, 
the quasiparticle energy and scattering rate of bosons 
have the following systematic expansions in the large- 
momentum limit (see Sec. II VI for details) 



E b (k) = 
and 



r>b 



1+64tt— f-r-2Re/ 
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-o(fc" 5 ) 



2m 
(2.4) 



r b (fc) 



0(k 



2m 
(2.5) 

Here, a\, is an s-wave scattering length between two iden- 
tical bosons, rib is a number density of bosons, and Cb 
is a contact density which measures the probability of 
finding two bosons close to each other [70j . f(\k\/n*) 
with being the Efimov parameter is a universal log- 
periodic function determined in Sec. IIVI [see Fig. [7] and 
Eq. (|4.31|) ]. Note that the coefficient of the contact den- 
sity in the scattering rate is always negative because Im/ 
ranges from —13.3 to —11.2. This rather counterintu- 
itively means that the energetic boson can escape from 
the medium easier than we naively estimate from a bi- 
nary collision. 

Each term has a simple physical meaning. Besides 
the free particle kinetic energy in Eq. (|2.2|) or (|2.4|) . the 
leading term represents a contribution from a two-body 
scattering in which the energetic atom collides with one 
atom coming from the atomic gas. The probability of 
finding such an atom in the atomic gas is quantified by 
the number density n a ^- Similarly, the subleading term 
represents a contribution from a three-body scattering in 
which the energetic atom collides with a small pair of two 
atoms coming from the atomic gas. The probability of 
finding such a small pair in the atomic gas is quantified 
by the contact density Cf j,. 

These results are valid for an arbitrary many-body 
state with translational and rotational symmetries (i.e., 
for any scattering length, density, and temperature) as 
long as Eq. (12.11) is satisfied. All nontrivial information 
about the many-body state is encoded into the various 
densities 7i CT .b and C/,6- 



B. Differential scattering rate 

We then consider the proposed scattering experiment 
in which we shoot a probe atom into the atomic gas and 
measure its differential scattering rate (see Fig. [T]). We 
assume that the probe atom denoted by x 1S distinguish- 
able from the rest of the atoms constituting the target 
atomic gas but has the same mass m, which is possible 
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by using a different atomic spin state. When the x atom ing rate has the following systematic expansion in the 
is shot into a spin-1/2 Fermi gas, its differential scatter- large-momentum limit (see Sec. [V] for details): 



^r x (fc) 



32 cos 9@ (cos 9) 



n t (x) + n±{x) 



Ifel 



32 1 2 cos 9 9(cos 9)k- <5(cos 9)k + 9(cos 9) p| 



Ifel 4 



2.9 cos 0, — 
v k' 



k\\ C f (x) 



0{k- 



*1 

2m 



(2.6) 



Here, 0( • ) is the Heaviside step function and 9 is a polar 
angle of the measured momentum p with respect to the 
incident momentum chosen to be k = \k\z. Accordingly, 
when a bunch of independent x atoms with a total num- 
ber N x is shot into the atomic gas, the number of scat- 
tered x atoms measured at an angle (9, ip) is predicted 
to be 

XJ.; V )-S X ^I*&&, (2.7, 

where the line integral is taken along a classical trajectory 
of the x atom. 

The differential scattering rate of the x atom shot into 
a spinless Bose gas is obtained from Eq. (12 .6[) by replacing 
the number density n^ + n^, the current density j'-j- + ji, 
and the contact density C/ of a spin-1/2 Fermi gas with 
those of a spinless Bose gas, rib, jb, and C&/2, respectively 
[see Eq. (|5.45p ]. These parameters are the same as those 
in Eqs. (|2.2[) - ([2.5[) , while translational or rotational sym- 
metries are not assumed here and thus the current density 
2a,b can be nonzero. On the other hand, n'^ in Eq. (12. 6[) is 
the Efimov parameter associated with a three-body sys- 
tem of the x atom with spin-up and -down fermions (the 
X atom with two identical bosons) and thus different from 
k* in Eqs. (|2.4p and (12.5[) associated with three identi- 
cal bosons. We note that the dependence on scattering 
lengths between the x atom and an atom constituting 
the target atomic gas appears only from 0(k~ 5 ) in the 
brackets. The corresponding formula in the weak-probe 
limit can be found in Eq. (I6.17[) . 

The first two terms in Eq. (|2.6p come from two-body 
scatterings and are proportional to the number density 
and current density of the target atomic gas. An impor- 
tant observation is that, because of kinematic constraints 
in the two-body scattering, they contribute to the for- 
ward scattering (cos 9 > 0) only. On the other hand, the 
last term comes from a three-body scattering and is pro- 
portional to the contact density. Its angle distribution is 
determined by a universal function g(cos 9, |fc|//«*), which 
is mostly negative on the forward-scattering side (see 
Fig. [9] in Sec. EJ. This is no cause for alarm, of course, 
because it is the subleading correction suppressed by a 
power of 1/| fe | to the leading positive contribution of the 
number density. 



In contrast, g(cos9, |fc|//c!|,) is positive everywhere on 
the backward-scattering side (cos 9 < 0) because it is 
now kinematically allowed in the three-body scattering. 
Therefore, the backward scattering is dominated by the 
contact density of the target atomic gas and its mea- 
surement can be used to extract the contact density in- 
tegrated along a classical trajectory of the probe atom 
[see Eq. (|2.7[) ]. Since the contact density is an impor- 
tant quantity to characterize strongly interacting atomic 
gases, a number of ultracold-atom experiments have been 
performed so far to measure its value but integrated over 
the whole volume 0, EH, El [HI Izl • Our proposed ex- 
periment provides a new way to locally probe the many- 
body aspect of strongly interacting atomic gases. 

Also we find from Eq. (|2.6p that the differential scatter- 
ing rate can depend on the azimuthal angle ip only by the 
current density of the target atomic gas. Therefore, the 
azimuthal anisotropy in the differential scattering rate 
may be useful to reveal many-body phases accompanied 
by currents. 



C. Comparison with Monte Carlo simulations 

All the above results are valid for an arbitrary many- 
body state at a sufficiently large momentum |fe| satisfy- 
ing Eq. (|2.ip . But how large should it be? One can gain 
insight into this question by comparing our results with 
other reliable results; for example, from Monte Carlo sim- 
ulations. Currently the only available Monte Carlo result 
comparable with ours is about the quasiparticle energy 
in a spin-1/2 Fermi gas [731 |74T|. Quite surprisingly, we 
will find reasonable agreement of our result (|2.2p with 
the recent quantum Monte Carlo simulation even at a 
relatively small momentum |fc|/fcp > 1.5. This indicates 
that our large-momentum expansions in Eqs. (|2.2p - (|2.6p 
are valid in a momentum range wider than we naively 
expect. 

In Ref . [ill , P. Magierski et al. extracted the quasipar- 
ticle energy E(k) from quantum Monte Carlo data in a 
balanced Fermi gas n-j- = nj_ at finite temperature. Their 
results are shown by points in Fig. [2] for (akp)^ 1 = 
at T/ep = 0.15 (left panel) and for (akp)- 1 = 0.2 at 
T/ep = 0.19 (right panel) in units of the Fermi energy 
ep = fc F /(2m) as functions of (\k\/kp) 2 . 
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FIG. 2. (Color online) Quasiparticle energies E(k)/tY as functions of (|fc|/fcp) 2 for (afcp) -1 = at T/e-p = 0.15 (left panel) and 
for (afcp)" 1 = 0.2 at T/ep = 0.19 (right panel). Points are results extracted from the quantum Monte Carlo simulation [7||] 
and dotted curves behind them are fits by a BCS-type formula (|2.13f> . Solid curves are our results from the large-momentum 
expansion (|2.8p with the use of the contact densities obtained in Refs. [75l [7r3 | . Narrow shaded regions behind them correspond 
to the contact densities varied by ±20% and broad ones indicate quasiparticle widths r(fc)/ep from the large-momentum 



expansion (12. 12[) with the same inputs. For comparison, free particle dispersion relations, Ef ICC (k) 
by dashed lines. 



fc7(2m) — /i, are shown 



Our large-momentum expansion of the quasiparticle 
energy (|2.2[) in the same units becomes 



E(k) 



jfe| 

ep \ kp 

+ o(k- 4 



16 1 ,r, C 

- rr- 7.54 -j 

o7T akp Kp 



(2.8) 



Here we used the definition of the Fermi momentum 
n a = fc F /(6"7r 2 ) and introduced the chemical potential 
\x because the quasiparticle energies in Ref. (74j are mea- 
sured with respect to /i. The values of chemical potential 
obtained in Ref. uM are 



— fa 0.471 and 0.319 



(2.9) 



for (akp)- 1 = at T/ep = 0.15 and for (akp)' 1 = 0.2 at 
T/ep = 0.19, respectively, with estimated errors of about 
10%. Note that the self-energy correction at 0[(fcp/|fc|) 2 ] 
receives two contributions: One is from the two-body 
scattering which can be attractive or repulsive depending 
on the sign of (akp) -1 . The other is from the three-body 
scattering which is proportional to C/kp > and thus 
always attractive. 

In order to make a comparison between our result 
and the Monte Carlo simulation, an input into the 
dimensionless contact density C/k F is needed ideally 
at the same scattering length and temperature as in 
Ref. [74|. The contact density at infinite scattering 
length (akp)' 1 = has been measured by a number of 
Monte Carlo simulations I7a-f79l | and ultracold-atom ex- 
periments [Tol [Til [Tgl [Til. |72|]. which are summarized in 
Table HI At zero temperature, they fall within the range 
of C/kp = 0.10 ~ 0.12. The temperature dependence 
of the contact density was reported in Refs. [Til [7§|, al- 
though the situation is somewhat controversial: The sim- 
ulation observed that the contact density increases with 



T/ep up to T/ep ps 0.4 while the experiment ob- 
served that the contact density monotonically decreases 
over the temperature range T/ep = 0.1 ~ 1 [Tij ]. Since 
the precise value of the contact density at T/ep = 0.15 is 
not yet available, we choose to use the value of Ref. [75[ 
at T/ep = 0.173(6); C/kp = 0.1102(11). This input fixes 
the self-energy correction to be 
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^37T akp k F J \ \k\ 

(2.10) 

and our result from the large-momentum expansion (|2.8[) 
is shown by the solid curve in Fig. [2] (left panel). In 
order to incorporate uncertainties of the contact density, 
its value is varied by ±20% which is represented by the 
narrow shaded region in the same plot. This variation of 
±20% is a very conservative estimate of the uncertainties 
because the contact density increases only by 15% even 
from T/ep ~ to 0.4 according to Ref. |79i |. 

On the other hand, the contact density away from the 
infinite scattering length is less understood, in particular, 
at finite temperature. Therefore, in order to facilitate 



TABLE I. Dimensionless contact density C/kp at infinite scat- 
tering length (afcp) - = from Monte Carlo simulations and 
ultracold-atom experiments at low and finite temperatures. 



Simulations 



Experiments 



Ref. 



C/fcp 



T/6F 



Ref. 



c/fcl r/6 F 



[77, 78] 
[76]" 

[79] 
[75] 
[79] 



0.115 

0.1147(3) 

0.0996(34) 

0.1102(11) 

0.1040(17) 







0.173(6) 
0.178 



[18] 0.118(6) 0.03(3) 
[10] 0.101(4) 0.10(2) 
[11] 0.105(8) 0.09(3) 
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a comparison between our result and the Monte Carlo 
result for (afcp) -1 = 0.2 at T/e F = 0.19, we choose to use 
the contact density of Ref. [76| for the same scattering 
length but at zero temperature; C/fcp = 0.156(2). This 
input fixes the self-energy correction to be 
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3-7T akF fcp J \ \k\ 

(2.11) 

and our result from the large-momentum expansion (|2.8p 
is shown by the solid curve in Fig. [5] (right panel). Note 
that the self-energy correction at (afcp) -1 = 0-2 is close 
to that at (a/cF)^ 1 = because opposite changes in the 
contributions from two-body and three-body scatterings 
happen to cancel each other. Again uncertainties of the 
contact density are incorporated by varying its value by 
±20% which is represented by the narrow shaded region 
in the same plot. 

In both cases of (afcp) -1 = and 0.2, one can see 
from Fig. [5] that our results are not very sensitive to the 
variations of the contact densities and, furthermore, they 
are in reasonable agreement with the quantum Monte 
Carlo simulation even at a relatively small momentum; 
(|fc|//cp) 2 > 2. This indicates that our large- momentum 
expansions in Eqs. (|2.2I) - (|2.6I) are valid in a wide range 
of momentum. 

Having our large-momentum expansions tested on the 
quasiparticle energy, we now present the quasiparticle 
width in the balanced spin-1/2 Fermi gas. Its large- 



momentum expansion 
becomes 



in units of the Fermi energy 



r(fc) 



16 kp _ 
' 37r]fe| ~ 
0{k'% 



1 

etfcp 



64 1 C 
37r ak-p kp 



(2.12) 



which is shown by the broad shaded region in Fig. [2] by 
using the same input into the contact density. Note that 
the correction at 0[(fcp /l^l) 3 ] vanishes for (o^f) -1 = 
(left panel), while it is given by -f-1.11 (fcjr/|fe|) 3 for 
(akp) -1 = 0.2 (right panel). In both cases, the quasipar- 
ticle widths gradually increase with decreasing momen- 
tum and eventually become comparable to the quasipar- 
ticle energies. This takes place at (|fe|/fcp) 2 ~ 2.1 and 
2.2, respectively, which roughly correspond to the point 
where our large-momentum expansions break down. 

P. Magierski et al. also extracted the pairing gap or 
pseudogap A, self-energy U, and effective mass m* pa- 
rameters by fitting a BCS-type formula 



^BCs(fc) 



V fc 2 



V 2m* 



■ U 



A 2 



(2.13) 



to their quasiparticle energies 74]. However, the fitted 
results (dotted curves in Fig. [2]) do not capture the cor- 
rect asymptotic behaviors at (\k\/kp) 2 > 3. This is be- 
cause i?BCs(fc) has the asymptotic expansion 

6f m* V Kf / £F £F 



in which the self-energy U < is taken to be a constant, 
while according to Eq. (|2.8[) , U should be momentum 
dependent and decay as 



► U — r - 7 - 54 JA 



ifei 



(2.15) 



at | k | / fep 3> 1 . Further analysis of their quantum Monte 
Carlo data incorporating our exact large- momentum ex- 
pansions may allow us better access to the intriguing 
pseudogap physics. 



III. SPIN-1/2 FERMI GAS 

Here we study properties of an energetic atom in a 
spin-1/2 Fermi gas and derive its quasiparticle energy 
and scattering rate presented in Eqs. (|2.2[) and (|2.3[) . 



A. Formulation 

The Lagrangian density describing spin-1/2 fermions 
with a zero-range interaction is 



C F = J2 V4 ( idt 



2m„ 



^+ci)\ip\i)^. (3.1) 



It is more convenient to introduce an auxiliary dimcr field 
4> = c ip^ip^ to decouple the interaction term: 



(3.2) 

For simplicity, we shall mainly consider the case of equal 
masses m = = m±. Some results in the case of un- 
equal masses are presented in Appendices [A] and |B] The 
propagator of fermion field in the vacuum is given by 



G(k) 



ko - £fe + i0 + 



2m 



(3.3) 



Also by using the standard regularization procedure to 
relate the bare coupling c to the scattering length a, 



If dk m m 
c /fcKA (2tt) 3 fc 2 47ra' 



(3.4) 



the propagator of dimer field cp in the vacuum is found 
to be 



D(k) 



-in 



1 



m k 2 i -n+ 1 

VT" mfc " - l0 a 



(3.5) 



D(k) coincides with the two-body scattering amplitude 
A(k) between spin-up and -down fermions up to a minus 
sign; A(k) = -D(k) (see Fig. [3). 
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v (i=i) 



FIG. 3. Two-body scattering amplitude iA(k) between spin- 
up and -down fermions. Solid (dashed) lines represent the 
fermion (dimer) propagator iG(k) [iD(k)]. Each fermion- 
dimer vertex (dot) carries i, and thus, iA(k) = — ID(fc). 



Our task here is to understand the behavior of the 
single-particle Green's function of spin-a fermions 



dye tky (T[^(x- 



!)]> 



(3.6) 



in the large-energy-momentum limit k — > oo for an arbi- 
trary few-body or many-body state. Without losing gen- 
erality, we can consider that of spin-up (a — f) fermions. 
The result for spin-down (a = \) fermions is obtained 
simply by exchanging spin indices f -l- 



B. Operator product expansion 



Im + 2.5 




25 mi 



FIG. 4. (Color online) Scaling dimension of the lowest three- 
body operator O3 = (m-j- + m_i)(j>(diipf) — mf-(di4>)ipf as a 
function of the mass ratio m^/mj_ taken from Ref. [831 . The 
solid curve is its real part and the dashed curve is its imagi- 
nary part shifted by +2.5. The scaling dimension of O3O3 is 
given by 2Re[Ae> 3 ]. 



for A 



o = % 

■^didjipa, -dityXdjipa), -didjtylipv), (3.12a) 



According to the operator product expansion [65 
l70l . H^-fUl ! the product of operators in Eq. (|3.6[) can be 
expressed in terms of a series of local operators O: 



dy e lky T[^(x + §)V4(z - 



!)] = 



E 



Wo t {k)Oi(x). 



(3.7) 

Wilson coefficients Wo depend on k — (k ,k) and the 
scattering length a. When the scaling dimension of a 
local operator is Aq, dimensional analysis implies that 
its Wilson coefficient should have a form 



W (k) = 



1 



fc|A +2 



Wo 



1 

oifci 



• ( 3 - 8 ) 



Therefore, the large-energy-momentum limit of the 
single-particle Green's function is determined by Wilson 
coefficients of local operators with low scaling dimensions 
and their expectation values with respect to the given 
state. 

The local operators appearing in the right-hand side of 
Eq. (|3.7[) must have a particle number No = 0. By recall- 
ing A^ ct —3/2 and A^ = 2 [87| . we can find twelve types 
of local operators with No = up to scaling dimensions 
A„ = 5: 



for A = 0, 



1 (identity) 



for A = 3, 



(3.9) 
(3.10) 
(3.11) 



iipld t ip a , id t (ipiipa), -i$di(t>, -id, 

(3.12b) 

for Ac = 5. Time-space arguments of operators (x) — 
(t, x) are suppressed here and below. Operators ac- 
companied by more spatial or temporal derivatives have 
higher scaling dimensions. 

Here we comment on scaling dimensions of operators 
involving more ip or <f> ~ ip^ip-x fields. Scaling dimensions 
of operators with three ip fields can be computed exactly 
by solving three-body problems (87l . |88| . For example, 
the lowest two operators are O = 2<p>(diip tJ ) — (di4>)i[i a 
and (jyipu and the products O^O have scaling dimen- 
sions A = 8.54545 and 9.33244, respectively. If more 
ip fields are involved, it is in general difficult to com- 
pute their scaling dimensions. However, with the help 
of the operator-state correspondence [U [H, |90| . they 
can be inferred from numerical calculations of energies 
of particles in a ha rmonic potential at infinite scatter- 
ing length 9ll4l00j |. For example, the ground-state en- 
ergy of four fermions, E — 5.01 hui [9(|, implies that the 
operator {<f>4>y (<f>4>) , which involves the lowest four-body 
operator <jxf>, has the scaling dimension A = 10.02. Be- 
cause adding more derivatives or fields generally increases 
scaling dimensions, we conclude that the operators in 
Eqs. (|3.9p ~ (|3.12p are the complete set of local operators 
with No = and Ac < 5 in the case of equal masses. 

In the case of unequal masses, however, this i s no t 
always the case due to the Efimov effect [lOll Il02j |. 
The scaling dimension of the lowest three-body oper- 
ator 03 = (m-j- + mj_)0(9iV't) — TUf(di(p)ipf decreases 
with increasing mass ratio / and eventually reaches 
Ao 3 = 5/2 at mf/m± = 13.607 so that A^tc, = 5 
(see Fig. [J). Furthermore, Ac 3 develops an imaginary 
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part for m^/m^ > 13.607, which indicates the Efimov ef- 
fect [H[ . In general, the Efimov effect for N particles im- 
plies that the corresponding TV-body operator On has the 
scaling dimension Aq n = 5/2 + is£ and thus A t q n = 5 
with si being a real number. The recent find ing of the 
four-body Efimov effect for rrif/mi > 13.384 [l02| indi- 
cates the existence of a four-body operator O4 whose scal- 



ing dimension becomes A 



\ 0i = 5 for m^jm\ > 13.384. 
Therefore, only when the mass ratio is below the low- 
est critical value for the Efimov effect, the operators in 
Eqs. (|3.9p - (|3.12[) are supposed to be the complete set of 
local operators with Nq = and A© < 5. 



C. Wilson coefficients 

The Wilson coefficients of local operators can be ob- 
tained by matching the matrix elements of both sides of 
Eg . (13.71) w ith respect to appropriate few-body states [H- 
70. |80h86| |. Details of such calculations are presented in 
Appendix [AJ In short, we use states (ip a {p')\ and \ip a (p)) 
to determine the Wilson coefficients of operators of type 
ipt4><j- The results are 



W t (k) = iG(k), 



(k) = -iG(k) 2 A(k), 



W 



1++ 



w 



(k) = - lG (k) 2 ^-A(k), 



1 d 2 



W 



t 



(*0 = ^G(k) 2 ^-A(k), 



dk 

G(k) 3 
4m 



(3.13) 
(3.14) 

(3.15) 
(3.16) 
(3.17) 



G(k) 



(3.18) 
(3.19) 

and all Wo(k) — for a =j, where A(k) is the two- 
body scattering amplitude between spin-up and -down 
fermions [see Eq. (|3.5[) ]. 

On the other hand, states (<t){p')\ and \4>(p)) are used 
to determine the Wilson coefficients of operators of type 
The results are 



= ^G(fc) 2 T t (£;,0;fc,0) 



(k) 



W 



dq 

(2tt) 3 

(fc)^ 



dq 



(3.20) 



' 2m J 



dq 




FIG. 5. Three-body scattering amplitude iT^(k,p;k',p') be- 
tween a spin- up fermion (solid line) and a dimer (dashed line). 



W_ 4# t^(*) = -iG(*) a ^r t (fc,p;*,p)| i ^ 



dq 



(3.21) 



dpi 



E P-b- 

2 ' 2 ' 



(3.22) 

Here, T^(k,p;k' ,p') is the three-body scattering ampli- 
tude between a spin- up fermion and a dimer with (k,p) 
\{k' ,p')\ being their initial (final) energy-momentum (see 
Fig. EJ. As we will show in Sec. UTTFl T t (k, 0; k, 0) and 
dT^(k,p; k,p)/dpi\ p ^o contain infrared divergences that 
are canceled exactly by the second terms in Eqs. (|3.20[) 
and (|3.2ip . respectively. Therefore, it is convenient to 
combine them and define finite quantities by 

T^(k, 0; k, 0) = T t (fc, 0; k, 0) - A(k) J A 



and 



(3.23) 



— T; ce (k,p;k,p)\ p ^ 



(3.24) 



Ar t (fc,p; k,p)\ p ^ - \JL A {k)j^L. (™ 



The regularized three-body scattering amplitude 
T t reg (fc,0;fc,0) will be computed in Sec. ImFl On the 
other hand, the ultraviolet divergences in the last 
two terms of Eq. (|3.20l) are canceled by those from 
expectation values of local operators as we will see 
below. 



D. Expectation values of local operators 

Now the single-particle Green's function of spin-up 
fermions for an arbitrary few-body or many-body state 
is obtained by taking the expectation value of Eq. (|3.7[) : 



dye ik y(T[iP t (x + §)V4(x - §)]) = ]T WoM(Oi(x)). 

(3.25) 
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The expectation values of local operators in Eqs. (|3.9 
(13. 12)) have simple physical meanings. For example, 



and 



(iptiJa) = n a (x) 



are the number density and current density of spin-u 
fermions and (^(V^Vv)) = din a (x), (didj^tp^)) = 

d t d 3 n a (x), (<9 t (V4Vv)) = d t n a (x), (#i(-iV£W<0) = 
dija(x) are their spatial or t emp oral derivatives. Fur- 
thermore, it is well known [68H701 |8fl-[Sa| that the expec- 
tation value of <f^<f) is related to the contact density C{x) 
by 



C(x) 



(3.28) 



function of spin-up fermions (|3.25[) is now written as 



(3.26) iG t (k)= dy^{T[^{x + ^\{x-l)]) 



(3.27) 



w t (k) + (k) ni + w_.^ (k) ■ j, 



+ w 



(k) 



+ W.. t M , (k) 



dq 

dq q 2 

(2tt) 3 2m 



r \ 1 c 



m- 



(3.32) 



By using the expressions of Wo[k) obtained in 
Eqs. (l3~T3"l) - (j3~2"2"j) . we find that Q t (k) can be brought 
into the usual form 



and {di{(j^(j))) = diC(x)/m 2 is its spatial derivative. 
The contact density measures the probability of finding 
spin- up and -down fermions close to each other [67H69] . 

jcj,{x) = m 2 {— i^'Vcf)) is an analog of the current density 
for dimer field 4> and shall be called a contact current 
density. 

If the given state is translationally invariant, the ex- 
pectation values of —ip^didjipa and iip' a dtip a can be ex- 
pressed in terms of the momentum distribution function 
of spin-cr fermions: 

pAl)= Jdye- i "-y{^ cr (t,x-^ (T (t,x+^)). (3.29) 

By using this definition, the expectation value of 
—ipldidjipa is found to be 



-ipldidji/) a } 



dq 



■<k<HPo{q)- 



(3.30) 



Similarly, the expectation value of iip^d^a can be eval- 
uated by using the equation of motion for fermion field 



(iipldtipa) 



(3 - 31) 



Both (—Tpldidji/ja) and (ixp^dtipo) contain ultraviolet 
divergences, which cancel those that already appeared in 
the last two terms of Eq. (pO0|) . 



E. Single-particle Green's function 

Since the derivatives of n a , j a , and C vanish for the 
translationally invariant state, the single-particle Green's 



g t (k) 



i 



k o -e k -^{k)+i0+' 



(3.33) 



where £-]-(&) is the self-energy of spin- up fermions given 
by 



d 



E t (k) = -A(k) H -—A(k)-j i 
1 92 A(k) t dq 



2 dkidkj 



(27T) 



3 



dq 



l^ A(k) 6 Ji 
2 dhdkj u 3 j (2tt) 3 



3 l ?<«7 

'" Q 2 [r (V)- 



dq q 2 ( C 



c 



Airrna 



^(fc,0;fc,0)^ - ^TT(*.K*.P)U ' 



TO 

(3.34) 



Here we eliminated the bare coupling c by using its re- 
lationship with the scattering length a [see Eq. (|3.4[) ]. 
By recalling the large-momentum tail of the momentum 
distribution function I67H69B 



lim p a {q) 

\q\-+oo q 



c 



0(q" 6 ), 



(3.35) 



one can see that the ultraviolet divergences in Eq. (|3.20p 
and Eqs. (|3.30|) and (|3.31|) canceled out so that S-)-(fc) is 
now manifestly finite. Corrections to the above expres- 
sion of S t ( fc ) denoted by "• • • " start with - (O) /k A °- 2 , 
where O are all possible operators with the lowest scal- 
ing dimension at Aq > 5. In the case of equal masses, 
A = 6. 

So far we only assumed that the given state is trans- 
lationally invariant. In addition, if the given state is 
rotationally invariant, j a , j^, and j dq/(2Tr) 3 (q i qj — 
Sijq 2 /3)p a (q) vanish. Therefore, in this case, the self- 
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energy of spin-up fermions is simplified to 

S t (fc) 



ir*(*,0iM)- ., 

TT1 ^ 



c_ 

m 2 



i—l 1 

dq q 2 

(2tt) 3 2^ 
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A(k) 



(3.36) 



Note that if the given state has the spin symmetry p^ = 
Pl, the integral of the momentum distribution function in 
the last line of Eq. p. 361) can be obtained from the energy 
density £ by using the energy relationship [67H69j | : 



E 



dq q 2 ( C 



= £- 



C 



Anma 



(3.37) 



The pole of the single-particle Green's function (13.33)) 
determines the quasiparticle energy and sc atter i ng r ate 
of spin- up fermions in a many-body system |l03l Il04| : 



k - e k - E-f(A;o,fc) = 0. 



(3.38) 



Because is as small as ~ 1/fc, we can set A;o = (-k 
in E-t-(fc) within the accuracy of (3(fc~ 4 ). Then the real 
part of the solution to Eq. (13.38)) gives the quasiparticle 
energy 

E f (k) = e k + Re[S t (e fe , fe)] + 0(fc~ 4 ), (3.39) 
while its imaginary part gives the scattering rate0 

r t (fe) = -2Im[£ t (e fe ,fe)] +0{k- i ). (3.40) 
Because 



dk ' 



0(fc- 4 ) (3.41) 



can be found by using the expression of A(k), we ar- 
rive at the following form of the on-shell self-energy of 
spin-up fermions for an arbitrary many-body state with 
translational and rotational symmetries^] 



S t (efc,fe) 



■in 



"I 



C 



i 'i 2 ma\k\ 3 

(3.42) 



4 0g (fc;fc)^ + 0(fc- 4 ). 



Here we denoted the regularized on-shell three- 
body scattering amplitude by t T * s (k;p) = 
T™ s (fc,0;p,fc-p)| fco = efc ,p o = £p . 

The first term in the on-shell self-energy (|3.42p is pro- 
portional to the two-body scattering amplitude A(ek, k) 
and the number density of spin-down fermions n\.. Its 
physical meaning is obvious: It is the contribution from 
the two-body scattering of the large-momentum spin-up 
fermion with a spin-down fermion in the medium. Simi- 
larly, the last term originates from the three-body scat- 
tering of the large- momentum spin-up fermion with a 
pair of spin-up and -down fermions close to each other, 
which is described by the dimer field cj>. The probabil- 
ity of finding such a small pair in the medium is given 
by the contact density C = m 2 (eft </>) [67l469lj. We note 
that the spin-down fermion and the small pair of spin- 
up and -down fermions coming from the medium are 
treated as being at rest because their characteristic mo- 
mentum ~ fep, A^ 1 are much smaller than |fe| in the large- 
momentum expansion [see Eq. (|2.1[) ]. 

Our remaining task is thus to determine the regular- 
ized on-shell three-body scattering amplitude i™ s (fc; fe) 
in Eq. (13.42)) up to 0(fc~ 4 ), which requires solving a 
three-body problem. Since it has a form of t™ s (k;k) = 

(m/fe 2 )if es [(a|fe|) -1 ], we need to determine the first two 
terms in its expansion in terms of (a, | fc | ) 1 . 



F. Three-body problem 



The quasiparticle residue Zf(k) is given by 



Z~\k) 



1 + 



4-7T 



ll 2 alfcl 4 



,-i-S t (fe) 

ni '-0(C/fc 4 ) 



Since Zf(k) = 1 + 0(k 4 ), the single-particle spectral density 
function of spin-up fermions becomes 



A t (k) 



-Im[S t (fc)] 



r t (fc) 



2n [ko - S f (fc)] 2 + [ir t (fc)]S 



within the accuracy we are currently working. Note that Af(k) 
at a large momentum \k\ ^ kp but below Fermi sea ko ~ — e^ (as 
opposed to ko ~ ej. in this paper) was studied in Ref . |105|| . The 
single-particle spectral density function was also computed in a 
self-consistent T-matrix approximation Il06ll and in a quantum 
cluster expansion at high temperature |107H . 



We now compute the three-body scattering ampli- 
tude Tf(k, 0; k, 0). Because Tf(k, 0; k, 0) does not solve 
a closed integral equation, we need to first consider 
T-|-(fc, 0;p, k—p), which is a solution to the integral equa- 
tion depicted in Fig. [51 and then take p — k. By denoting 



3 In the case of unequal masses m-f ^ m^, this result is modified 
to 

4tt nj i m-fC 



E+(e&f ,fc) 



r H_ 

~~M ' aTfcT 



m i , _!_ 

hi T a\k\ 



(2^) 2 a|fe| 3 



-T5 eg (fc,0;fc,0)|, -±_ + 0(fc- 4 ), 

where = fe 2 /(2m„), a total mass M = m-j- + mj_, and a 
reduced mass fj, = m^m^/(m^ + m±). 
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FIG. 6. Integral equation for the three-body scattering amplitude T^(k;p) = T^(k,0;p,k—p) between a spin-up fermion and a 
dimer. 



T^(k,0;p,k—p) simply by T^{k;p), its integral equation 
is written as 



T t (A;p) - G(-p) 
*dq Q dq 



T^k;q)G{q)D{k-p)G{k-p-q). 



(3.43) 



Because Tf(/c; q) is regular in the lower half plane of qo, 
the integration over qg can be easily performed to lead to 



T f (k;p) = G(-p) 
dq 



(20 



3 T t (k;q)D(k-q)G(k-p-q)\ 



(3.44) 



qa=tq 



Then by setting fco = po = £p an d defining tf(k;p) = 
Tf(ek, fe; £p,p), we obtain an integral equation solved by 
the on-shell three-body scattering amplitude 

771 f dq 

t t (k;p) = — V rlC a (k;p,q)t t (k;q), (3.45) 

where the integral kernel K a {k\p, q) is defined by 

4tt 



1 



(3.46) 



q 2 +pq kp kq iCH 



This integral equation has to be solved numerically 
to determine t^(k;p). We have computed t^(k;p) at 
(a|fe|) _1 = and its first derivative with respect to 
(a. | fc | ) ~ 1 . More details of solving the integral equation 
(I3.45[) are presented in Appendix [5] 

As we mentioned after Eq. (|3.22p . tf(k;k) — 
T^(k, 0; k, 0)\k o =e k contains an infrared divergence. This 
can be seen by rewriting p.45[) at p = k as 



t t (k;k) 



/C a (fc;fc,g)_ 



(27T) 



dq 



K-a(k;k,q) 



t t (k;q) + -z 



(3.47) 



in which the second term is infrared divergent because 

4tt 1 



K. a (k;k,q) 



I q 2 



(3.48) 



at \q\ — > 0. However, this infrared divergence is can- 
celed exactly by the second term in Eq. (|3.23p . Therefore, 



the regularized three-body scattering amplitude defined 
there is free of divergences and its on-shell version is given 
by 



4 cg (fc;fc)=T t rog (M;fc,0)k =efc 



rn 



dq 



JC a (k;k,q) 



47T 



;|fe| 



1 q 2 



rn 
1? 



dq 



JC a (k;k,q) 



Tfl 



(3.49) 



By using the numerical solutions of tf(k; q) at (<x|fc|) 1 = 
and its first derivative, the expansion of t™ e (k;k) in 
terms of (a|fc|) _1 is found to be 



3.771 + l^z + 0(fc- 2 ) 
are 



m 



(3.50) 



These numbers are universal (i.e., independent of short- 
range physics). We note that the Born approxima- 
tion [the first term in Eq. (|3.49|) ] gives t™ 6 (k; fe)|Bom = 
(—1) m/fc 2 , which is wrong even in its sign. 

The necessity of the above subtraction procedure is 
physically understood in the following way: The "bare" 
three-body scattering amplitude tf(k;k) describes the 
three-body scattering of the large-momentum spin-up 
fermion with a pair of spin-up and -down fermions at rest. 
tf(k;k) is obtained from the integral equation depicted 
in Fig. [6j which actually includes a process in which the 
large-momentum spin-up fermion collides only with the 
spin-down fermion coming from the pair and the other 
spin-up fermion remains a spectator staying away from 
the scattering event. This is seen by recalling that the 
momentum q in the second term of Eq. ()3.47j) is the rel- 
ative momentum between spin-up and -down fermions 
constituting the pair and thus small \q\ corresponds to 
the large interparticle separation. This process is essen- 
tially the two-body scattering which is already included 
in the first term of Eq. (|3.42[) . Therefore, to avoid the 
double counting, the contribution of this two-body scat- 
tering process has to be subtracted from the bare three- 
body scattering amplitude. This leads to the regularized 
three-body scattering amplitude t™ s (k;k) in Eq. (|3.49p . 
which appears in front of the contact density in the last 
term of Eq. 

Finally, by substituting the numerical solution of the 
three-body problem (|3.50p into Eq. (|3.42p . we find that 
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the on-shell self-energy of spin-up fermions has the follow- 
ing systematic expansion in the large-momentum limit: 



16tt 



7.54 



2 

a\k 
22.1 . 
4fcf 



a 2 \k\ 2 
C 

W 



Ifel 3 



0(fc" 6 ) 



(3.51) 



This result combined with Eqs. ([339)1 and (|3~40| leads 
to the quasiparticle energy and scattering rate of spin-up 
fermions presented previously in Eqs. (|2.2I) and 



IV. SPINLESS BOSE GAS 

Here we study properties of an energetic atom in a 
spinless Bose gas and derive its quasiparticle energy and 
scattering rate presented in Eqs. (|2.4j) and (|2.5p . The 
analysis is similar to the previous case of a spin- 1/2 Fermi 
gas. 



A. Formulation 

The Lagrangian density describing spinless bosons 
with a zero-range interaction is 



C u = ^ [ idt + ^ + 1 
2m / 2 



(4.1) 



/ V 2 \ 1 1 1 

= ijy I idt + — J if) — —fcj) + -(piM + -^tyV> 
\ 2m ) Ac I 2 

where an auxiliary dimer field 4> = ctpip is introduced to 
decouple the interaction term. The propagator of boson 
field i\) in the vacuum is given by 



G{k) 



1 



ko - ek + i0 + 



fc_ 

2m 



(4.2) 



Also by using the standard regularization procedure to 
relate the bare coupling c to the scattering length a, 



1 



dk m m , 

~7~' ( 4 -3) 



fe|<A (2tt)3 fe2 47ra' 



the propagator of dimer field 4> in the vacuum is found 
to be 



D(k) 



Sn 



(4.4) 



D(k) coincides with the two-body scattering amplitude 
A(k) between two identical bosons up to a minus sign; 
A(k) = -D(k) (see Fig. |3J|. Note that D(k) and A{k) in 
the case of bosons are twice as large as those for fermions 
[see Eq. (1531) ]. 



Our task here is to understand the behavior of the 
single-particle Green's function of bosons 



dye iky m{x+m^x- §)]) (4.5) 



in the large-energy-momcntum limit k — > oo for an arbi- 
trary few-body or many-body state. 



B. Operator product expansion 

According to the operator product expansion [68 
70. 80 86]. the product of operators in Eq. (ITU can be 
expressed in terms of a series of local operators O: 

Jdy^ ky T[i>{x+^\x-fj\=Y,W 0i {k)Oi{x). 

(4.6) 

The local operators O appearing in the right-hand side 
must have a particle number Nq = 0. By recalling A^, = 
3/2 and = 2 [83], we can find thirteen types of local 
operators with No = up to scaling dimensions Aq = 5: 



for A = 0, 



for A = 3, 



1 (identity) (4.7) 
ipty (4.8) 
i^d^, -tfi,(^ty), <t>H (4.9) 



for A 



o 



■tfdidrf, -diitfdrf), -didj{^), (4.10a) 



i^dti/;, idti^H), -itfdi<t>, -t&C^V). (4.10b) 



(4.10c) 



for Aq — 5. 

The striking difference between the cases of fermions 
and bosons is the presence of the Efimov effect in a sys- 
tem of three identical bosons [l08f . As we discussed at 
the end of Sec. IIIIB1 the Efimov effect implies that the 
corresponding three-body operator tfnp has the scaling di- 
mension A = 5/2 + iso so that ((/rip )^(<^) in Eq. (|4.10c[) 
has A = 5 [lO, |88| . The determination of its Wilson co- 
efficient requires solving a four-body problem, which is 
beyond the scope of this paper. Therefore, in this sec- 
tion, we only consider the local operators with Nq = 
and A < 4 in Eqs. (|4.7j) - (l4.9j) . As before, their expec- 
tation values have simple physical meanings such as the 
number density of bosons, 



(ijj'ip) = n(x), 



(4.11) 
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and its spatial derivative (di^ tp)) 
density of bosons, 



difiix), the current 



(-i^Vip) = j(x), 



and the contact density 



C{x) 



(4.12) 



(4.13) 



which measures the probability of finding two bosons 
close to each other. This definition of the contact density 
coincides with that used in Ref. [70j and thus the large- 
momentum tail of the momentum distribution function 
of bosons is given by limiqi^^ p(q) = C/q 4 + 0(q~ 5 ). 



C. Wilson coefficients 

The Wilson coefficients of local operators can be ob- 
tained by matching the matrix elements of both sides of 
Eq . (I4.6|) w ith respect to appropriate few-body states (68l — 
ProL l80l - [86| . Details of such calculations are presented in 
Appendix [S] The results are formally equivalent to those 
in the case of fermions as long as A(k) in Eqs. f|3 . 13|) — 
p.22j) is understood as the two-body scattering ampli- 
tude between two identical bosons [see Eq. 



W t (k)=iG(k), 



W^(k) = -iG(k) 2 A(k), 



W_ i ^Jk) = -iG(kf-^A(k), 



W_ idz (^^)(k) = 0, 



W^4k) = -iG(k) 2 T(k,0;k,0) 



dq 

J2n) 3 \q 2 



W^t^(fc) 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



Here, T(k,p; k' ,p') is the three-body scattering ampli- 
tude between a boson and a dimer with {k,p) [(£;', p')] 
being their initial (final) energy- momentum (see Fig. [5J . 
Because T(fc, 0; fc, 0) contains an infrared divergence that 
is canceled exactly by the second term in Eq. (|4.18p . it is 
convenient to combine them and define a finite quantity 

by 



T reg (fc, 0; k, 0) = T(k, 0; k, 0) - A{k) 



dq 

(2^)" 3 \q 2 



(4-19) 

This regularized three-body scattering amplitude will be 
computed in Sec. IIVEI 



D. Single-particle Green's function 

Now the single-particle Green's function of bosons for 
an arbitrary few-body or many-body state is obtained by 
taking the expectation value of Eq. (|4.6|) . By using the 
expressions of Wo(k) obtained in Eqs. (|4.14[) - (14.18l) . we 
find that it can be brought into the usual form 



iG{k) = dye lk y(T[^(x+^(x - §)]) 



(4.20) 



k -e k -T,{k)+iQ+' 
where S(fc) is the self-energy of bosons given by 

E(fc) = -A(k) n - ^-A(k) ■ j - T rog (fc, 0; k, 0)-^ - • • • . 

ok ra z 

(4.21) 

Corrections to this expression denoted by "• • • " start 
with ~ (0}/k A °~ 2 , where O are all operators in 
Eq. (|4.10p with the scaling dimension Aq = 5. 

As in Eq. (|3.38|) , the pole of the single-particle Green's 
function (|4.20[) determines the quasiparticle energy and 
scattering rate of bosons in a many-body system. Within 
the accuracy of 0(k~ 4 ), the real part of£(efc,fc) gives the 
quasiparticle energy 

E{k) = e k + Re[S(e fc , fe)] + 0(fc~ 4 ), (4.22) 

while its imaginary part gives the scattering rate 

r(Jfe) = -2 Im[S(e fe , fc)] + 0(fc~ 4 ). (4.23) 



By setting fco = e& in Eq. (I4.21j) . the on-shell self-energy 
of bosons for an arbitrary state is found to be 



S(cfc,fe) 



8tt 



Airi 



k 



t -g (fc;fc) ±L_ 

Ul 



Q{k 



1 i 1 1 

2 a\k\) 



m\k\ 



(4.24) 



where we denoted the regularized on-shell three- 
body scattering amplitude by i rcg (fc;p) = 
T rcg (fc ,0-p , fc-p)|fc =e fc , Po ^e p . The second term in 
Eq. (|4.24l) . which is proportional to dA(k , k) /dk\k ^ Ch 
and the current density ji, represents the contribution 
from the two-body scattering in which the large- 
momentum boson collides with a boson moving with a 
small momentum. The physical meanings of the other 
two terms were discussed at the end of Sec. IIII El 

Our remaining task is thus to determine the regular- 
ized on-shell three-body scattering amplitude t ICS (k; k) 
in Eq. (|4.24p up to 0(fc~ 3 ), which requires solv- 
ing a three-body problem. Because three identical 
bosons suffer from the Efimov effect, (k 2 /m) i rcg (fc; k) = 
i rcg [(a|fc|) _1 , |fc|//c*] depends not only on (a|fc|) _1 but 
also on I/cI/k*, where is the Efimov parameter. As long 
as we are interested in t r ° s (k; k) ~ m/fe 2 within the accu- 
racy of 0(k~ 3 ), we can set the scattering length infinite 
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(a|fc|) _1 = 0, because the dependence on it appears only 
from 0(k~ 3 ). Then the resulting quantity f rog [0, |fc|//c»] 
needs to be determined, which is a log-periodic function 
of |fe|//«* as we will see below. 



E. Three-body problem 

We now compute the three-body scattering ampli- 
tude T(k, 0; k, 0). Because T(k,0;k,0) does not solve 
a closed integral equation, we need to first consider 
T(k,0;p,k—p), which is a solution to the integral equa- 
tion depicted in Fig.El and then take p = k. By denoting 
T(k,0;p,k—p) simply by T(k;p), its integral equation is 
written as 

T(k;p) = -G(-p) 

+ ij^T(k; q )G( q )D(k- P )G(k- P - q ). (425) 

Then by performing the integration over qo and defining 
t(k;p) = T(efc, k; e p , p), we obtain an integral equation 
solved by the on-shell three-body scattering amplitude 

m ^ 9 f dq 



t(k;p) 



K a {k;p,q)t{k;q), (4.26) 



p' J (2tt) 3 

where the integral kernel /C a (fc;p, q) is defined in 
Eq. (|3.46[) . This integral equation has to be solved nu- 
merically to determine t(k;p) at (a|fc|) _1 = 0. 

Compared to the case of fermions in Eq. (|3. 451) . both 
signs in the right-hand side of Eq. (|4.26[) are opposite due 
to different statistics and the second term has the fac- 
tor 2 originating from the fact that the two-body scat- 
tering amplitude between two identical bosons is twice 
larger. The former di fference l eads to the striking conse- 
quence: As is known |l09l llld |. the three-boson problem 
described by the integral equation (|4.26[) is ill defined 
without introducing an ultraviolet momentum cutoff A 
in the zero orbital angular momentum channel. A can be 
related to the Efimov parameter which is defined so 
that three identical bosons at infinite scattering length 
have the following infinite tower of binding energies: 



2 

-2-rrn/so 



[n 



oo 



(4.27) 



with s = 1.00624 [ll0llTTl| . Because K * IS defined up to 
multiplicative factors of A = e*l Sa = 22.6944, the solu- 
tion to the integral equation (|4.26[) at (a|fe|) _1 = has to 
be a log-periodic function of |fe|//t*. We have computed 
such t(k;p) in a range 1 < |fc|/«* < A 2 corresponding to 
two periods. More details of solving the integral equation 
(|4.26[) are presented in Appendix [B] 

As we mentioned after Eq. (|4.18[) , t(k;k) = 
T(k, 0; k, 0)|fc o=efc contains an infrared divergence. This 
can be seen by rewriting (|4.26[) at p = k as 

— —JC a {k;k,q) — 



dq 



(27T) 

K. a (k;k,q) 



t(k;q)- — 
9 



(4.28) 




FIG. 7. (Color online) Universal log-periodic func- 
tion /(|A;|/k») defined in Eq. (|4.30[) as a function of 
(so/i") ln(|fe|//t*). Circles (squares) with steps 1/16 corre- 
spond to its real part (imaginary part shifted by +12) and 
solid curves are fits by the approximate formula (|4.31f) . Two 
periods in the range 1 < |fe|/«* < e 2lT ^ S0 are shown here. 



in which the second term is infrared divergent at | q \ — > 
[see Eq. (|3.48p ]. However, this infrared divergence is can- 
celed exactly by the second term in Eq. (|4.19p . Therefore, 
the regularized three-body scattering amplitude defined 
there is free of divergences and its on-shell version is given 



by 



f^(k- 1 k)=T^(k,0;k,0)\ ko=eh 
dq 



m 

V +2 



(27T) 



K, a (k;k,q) 



-in 



dq 



K-a(k;k,q) 



t(k;q) 



2 a 



m 

7^2 



in 
T2 



9" 

(4.29) 



The physical meaning of this subtraction procedure was 
discussed at the end of Sec. IIII Fl 

By using the numerical solutions of t(k; q) at 
(a|fe|) _1 = 0, t TCg (k;k) is computed in the range 1 < 
|fe|//t* < A 2 and the resulting universal function 



/(^U-r*(fc;fc)| 
At* / m 



(a| fe|) — 1=0 



(4.30) 



is shown by points in Fig. [71 The logarithmic periodicity 
is clearly seen and we find that our numerical results are 
excellently reproduced by 



f(z)~X + 



Yi cos(2sq In z + 6i) + %Y% sin(2so lnz + 6\) 



1 + Zsm(2s lnz + 5 2 ) 

(4.31) 

with fitting parameters X « -0.09656 - 12.20 i, Y x w 
1.036, Y 2 w -1.032, Z « -0.08460, and 8 X w S 2 w 
0.4653. The approximate formula (14.31)) is plotted by 
solid curves in Fig. [7j and differs from the numerical 
points only by the amount < 4 x 10~ 6 . Whether 
Eq. (|4.31l) is the true analytic expression of /(|fc|/ft„) 
or not needs to be investigated further. 

Finally, by substituting the numerical solution of the 
three-body problem (|4.30p into Eq. (|4.24|) . we find that 
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the on-shell self-energy of bosons has the following sys- 
tematic expansion in the large- momentum limit: 



E(c fc ,fc) 



32tt -i 



a\k\ 
k\\ C 



- 32m 



k- j 



(4.32) 



<"fc- 



This result combined with Eqs. (I4.22|) and (|4.23|) leads 
to the quasiparticle energy and scattering rate of bosons 
presented previously in Eqs. (|2.4I) and (12.51) . where the 
contribution of the current density is dropped by assum- 
ing translational and rotational symmetries. 



DIFFERENTIAL SCATTERING RATE 



Bose gas with a large momentum and measure its angle 
distribution. The differential scattering rate presented in 
Eq. (I2.6[) will be derived from the total scattering rate by 
using the optical theorem, while it coincides with the one 
expected on physical grounds. Here translational or rota- 
tional symmetries are not assumed and thus the densities 
depend on a time-space coordinate (x) = (t,x). 



A. Spin- 1/2 Fermi gas 

We first consider the case of a spin- 1/2 Fermi gas. Here 
a probe atom is denoted by x and assumed to interact 
with spin-up and -down fcrmions by scattering lengths a-f- 
and oj,, respectively. The Lagrangian density describing 
such a problem is 



So far we have studied a quasiparticle energy and a 
"total" scattering rate of an energetic atom both in a 
spin-1/2 Fermi gas (Sec. Mil) and in a spinless Bose gas 
(Sec. ITVj) . Often in physics, differential scattering rates 
or cross sections also reveal many important phenom- 
ena. For example, differential cross sections in neutron- 
deuteron or proton-deuteron scatterings at intermediate 
or higher energies are important to reveal the existence 
of three- nucleon forces in nuclei [29M3T| . Also, momen- 
tum and angular resolutions have been essential to re- 
veal short-range pair correlations in nuclei from two- 
nucleon knockout reactions by high-energy protons or 
electrons [32|,[33j]. Furthermore, differential cross sections 
of high-energy neutrons scattered by liquid helium have 
been employed to extract the momentum distribution of 
helium atoms (lol - IZij] . 

Now in ultracold-atom experiments, one can in prin- 
ciple imagine shooting an energetic spin-up fermion (bo- 
son) into a Fermi (Bose) gas trapped with a finite depth 
and measure the angle distribution of spin-up fermions 
(bosons) coming out of the trap. Here one needs to 
be cautious, however, because the incident atom can- 
not be distinguished from atoms constituting the atomic 
gas. For example, there is a process in which the ener- 
getic spin-up fermion collides with a spin-down fermion 
in the medium and they escape from the trap. However, 
such a spin-down fermion may be accompanied by an- 
other spin-up fermion nearby so that they form a small 
pair described by the contact density. What happens to 
this spin-up fermion when its partner is kicked out by 
the incident atom? Whether it escapes from the trap 
to be measured or not has to be imposed consistently 
on all calculations, which appears intractable in our sys- 
tematic large-momentum expansion without introducing 
phenomenological procedures. 

In order to avoid this problem and unambiguously de- 
termine the differential scattering rate, it is therefore fa- 
vorable to consider an incident atom that is distinguish- 
able from the rest of the atoms constituting the atomic 
gas. In this section, we imagine shooting a different spin 
state of atoms into a spin-1/2 Fermi gas or a spinless 



2rn-i 



- £ ^xHt^aX + ^-F, (5.1) 

<7=t4 



where Lp defined in Eq. (|3.ip describes spin-1/2 fermions 
interacting with each other by a scattering length a. For 
simplicity, we shall assume that all particles have the 
same mass m = m x = = m±. In analogy with 
Eq. (13. 5p , the two-body scattering amplitude between the 
X atom and a spin-er fermion is given by 



A a (k) 



47T 



1 



(5.2) 



The behavior of the single-particle Green's function of 
the x atom, 

iG x (x;k) = Jdy^{T[x{x+^)x\x-l)]), (5.3) 

in the large-energy-momentum limit k — > oo can be un- 
derstood by using the operator product expansion as in 
Sees. |HT] and [IVJ Since three distinguishable particles 
(Xi t) i) wit h zero-ran ge interactions suffer from the Efi- 
mov effect jl lOl llll| . we only consider local operators 
with No = up to scaling dimensions Aa — 4 (see dis- 
cussions m Sec. ITVBj) . Then in analogy with Eqs. (TOU1) 
and (|4.21l) . G x (x; k) can be written in the usual form 



G x (x;k) = 



1 



h) - (k ~ £ x (a;; k) + iO H 



(5.4) 



where S x (a;; k) is the self-energy of the x atom given by 



S x (a;; k) 



<T=t4 



•3?*(M;M) 



Ar(fc) n a (x) 
C(x) 



d_ 

dk 



A a {k)-j a {x) 



(5.5) 



Here, n a (x) — (if)^if> a ) and j a {x) = (—iiplVijj a ) are the 
number density and current density of spin-tr fermions 
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FIG. 8. Set of integral equations for three-body scattering amplitudes T x (k;p), Tf(k;p), and Tj,(fc;p) involving three distin- 
guishable particles \ and spin-t and -\, fermions. Momentum labels are the same as those in Fig. [6] 



and C(x) — m (<j>*<j>) is the contact density of a spin- 
1/2 Fermi gas. Note that these parameters only refer 
to many-body properties of the given spin- 1/2 Fermi gas 
and do not involve the information related to the \ atom. 
On the other hand, T x s (k, 0; k, 0) is a finite quantity 
defined by 



I£*(fc,0;M) 

= T x (k,0;k,0)-[A t (k)+A^(k)] 



(5.6) 



dq 

W) 3 W 



where T x (k,p; k! ,p') is the three-body scattering ampli- 
tude between the x atom and a dimer composed of spin- 
up and -down fermions with (k,p) [(k',p')] being their 
initial (final) energy-momentum (see Fig. [5]). Corrections 
to the above expression of ^ x (x; k) denoted by "• • • " start 
with ~ (0}/k A °~ 2 , where O are all possible operators 
with the scaling dimension Aq = 5. 

Then by setting k = e k in Eq. (|5.5I) . the on-shell self- 
energy of the x atom for an arbitrary state is found to 
be 



S x (x;e fc ,fc) 



47r n a {x) 



2-n-i k-j a (x) 



a„\k\ 



m\k\ 



i ) 



m\k\ 



f^(k;k)^4 + 0(k- 3 ). 



(5.7) 



The quasiparticle energy and scattering rate of the x 
atom in a spin- 1/2 Fermi gas are given by the real and 
imaginary parts of ^ x (x; e k , k) according to 



E x (x; k) = e k + Re[E x (x; e k , k)} + 0(fc" 4 ) (5.8) 



and 



T x (x; k) = -2 Im[E x (x; e k , k)} + 0(fc -4 ), 



(5.9) 



respectively. Our next task is to determine the regular- 
ized on-shell three-body scattering amplitude t T x s (k; p) = 
T x eg (fc,0;p, fc-p)|fc = efc ,p =e p in Eq. (1B"J|) . which requires 
solving a three-body problem. 



B. Three-body problem 



We now compute the three-body scattering amplitude 
T x (k, 0; fc, 0). Unlike the previous cases in Sees. IIIIFI and 
ITTEl T x (k;p) = T x (k,0;p,k-p) by itself does not solve 
a closed integral equation. To find a closed set of integral 
equations, we need to introduce other three-body scatter- 
ing amplitudes T a (k;p) with a = f, J,, which describe pro- 
cesses where the x atom and a dimer composed of spin-up 
and -down fermions with their energy-momentum (fc,0) 
are scattered into a spin-<r fermion and a dimer com- 
posed of the x atom and the other fermion with their 
energy-momentum (p,k—p). These three scattering am- 
plitudes are solutions to a closed set of integral equations 
depicted in Fig. [5J Then by following the same pro- 
cedures as in Eqs. (|3.43[) ~ (13.45[) . the integral equations 
solved by the on-shell three-body scattering amplitudes 
t x ^^(k;p) = T x ^^(e k , k; e p ,p) can be written as 

*x( fe ;p) = J ^^£a 4 (fc;p,<z)if(fc;g) 



(27T) 

dq 



(5.10a) 



/C QT (fc;p, q) t±(k;q), 



h{k;p) = ^ + I ' j7T^fca(k;p,q)t x (k;q) 

p J {^r 



dq 

(2^y 



(5.10b) 



JC at (k;p, q) t±(k;q), 
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771 f d>Q 

h( k >P) = ~2 + / 77rT3^a(fe;p,q)t x (A;;q) 
p- J {Zir) 3 



dq 



(5.10c) 



JC ai {k;p, q) t^{k;q) 



where the integral kernel K, a (k;p,q) is denned in 
Eq. GOB)- 

As long as we are interested in t^ s (k;k) ~ m/k 2 up 
to 0(k~ 3 ) [see Eq. (|5.7p ], we can set all three scatter- 
ing lengths infinite (a-|-|fc|) _1 , (aj.|fc|) _1 , (a|fc|) _1 = 0, 
because the dependence on them appears only from 
0(k~ 3 ). In this case, by defining 



t F (k;p) 



and 



t B (k;p) 



2t x {k;p) - t t (k;p) - ti(k;p) 



t X (k;p) +t t (k;p) + t;(fc;p) 



(5.11) 



(5.12) 



the three coupled integral equations (|5. 10[) can be 
brought into two independent integral equations: 



t F (k;p) = j 

and 



t B (k;p) 



P 



dq 

(2^)3 

dq 



K, 00 (k;p,q)t F (k;q) (5.13a) 



(2tt)* 



K 00 {k;p,q)t B {k;q). (5.13b) 



Because these two integral equations are equivalent to 
Eqs. (f3T45|) and ({4726]) at (o,| A;| )^ 1 = 0, their solutions 
are already obtained. 

By using the numerical solutions obtained previously 
in Eqs. (13750)) and (t4730|i . £™ g (fc;fc) within the accuracy 
of 0(k~ 3 ) is found to be 



and 

T x (x;k) 




(5.16) 



respectively. 



C. Differential scattering rate 

Our final task is to determine the differential scattering 
rate dT x (k)/dp; that is, the rate at which the x atom 
shot into a spin- 1/2 Fermi gas with the initial momentum 
k is measured at the final momentum p (see Fig. [IJ. We 
start with the following expression for the total scattering 
rate obtained from Eqs. (I5.7[) and ([52 



T x (x;k) 

= E 



4tt 



n a (x) 



2Im^(fc;fe)^M 

<7=t,4. 



^ (q^lfcl) 2 k-j a (x) 
'1 | i I 2 mlfcl 2 



<3(fc 



-3\ 



4 n -)(jfe) + rW")(jfe) +r( c )(fc) + o(fc 



,-3\ 



(5.17) 

Here the contributions of the number density n CT! current 
density and contact density C are denoted by T 



r^'\ and respectively. 



,(C) 



Contribution of number density 



f x e %k; k) = - [t£*(fc; k) + t r ^(k; k)] + 0(fc~ 3 ) 



3.771 + / ^ 



£ + 0(k- 3 ). 



(5.14) 



Here, f(\k\/K' m ) is the universal log-periodic function 
plotted in Fig. [7] and approximately given by Eq. (|4.31[) 
with k'x being the Efimov parameter associated with a 
three-body system of the x atom with spin-up and - 
down fermions. By substituting this numerical solution 
of the three-body problem into the on-shell self-energy in 
Eq. (|5.7[) . the large- momentum expansions of the quasi- 
particle energy and scattering rate of the x atom in a 
spin-1/2 Fermi gas are obtained from Eqs. fl575} and 



E x (x;k) 



C(x) 
Ifcl 4 



+ 0{k- 5 ) 



(5.15) 



In order to extract the differential scattering rate, it is 
instructive to rewrite the first term in Eq. (|5.17p so that 
its physical meaning becomes transparent: 



[7"U) -».ri.r) f^\A a (e k ,k)\ 2 



(5.18) 



x (2n) i 5(p + q — k)S(e p + t q - e k ). 

Here, A a (k) introduced in Eq. (|5.2[) is the two-body 
scattering amplitude between the x atom and a spin- 
a fermion. It is now obvious that ^(fe) represents 
the contribution from the two-body scattering in which 
the x atom with the initial momentum k and a spin-u 
fermion at rest are scattered into those with their final 
momenta p and q, respectively. Therefore, the contri- 
bution of the number density of spin-cr fermions to the 
differential scattering rate of the x atom can be read off 
as 



dY ( x a \k) 
dp 



= n a (x) 



dq 



\A a (e k ,k) 



(2tt) 6 1 

x (2n) i 5(p + q- k)S(e p + e q 



(5.19) 



Efe) 
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Then by performing the integration over the magni- 
tude of momentum |p| 2 <i|p|, the angle distribution of the 
scattered \ atom is found to be 



dT { x "\k) _ 4cos0 9(cos0) n CT (a;) 



(HI 



\k\ 



(5.20) 



where 9( • ) is the Heaviside step function and is a polar 
angle of the final momentum p with respect to the initial 
momentum chosen to be fc = \k\z. Because of kinematic 
constraints (energy and momentum conservations) in the 
two-body scattering, the number density contributes to 
the forward scattering (cos > 0) only0 



2. Contribution of current density 



Similarly, the second term in Eq. (|5.17|) can be rewrit- 
ten as 



r ^ ,(fc) ^ ;w i/H iA ^ fc) 



(27r) 4 <5(p + q k)S(e p + e q - ko) 



(5.21) 



which represents the contribution from the two-body 
scattering in which the \ atom is scattered by a spin- 
cr fcrmion moving with a small momentum. Accordingly, 
the contribution of the current density of spin-cr fermions 
to the differential scattering rate of the x atom can be 
read off as 



dT 



(jV), 



k) 



dq 



dp Juy ~' dk J (2?r) 

x (2ir) 4 5(p + q — k)S(e p + e q - k Q 



6 \A a (k ,k)\ 2 



(5.22) 



ko— >£fc 



Then by performing the integration over the magni- 
tude of momentum |p| 2 <i|p|, the angle distribution of the 
scattered \ atom is found to be 

dT { 4°\k) _ 2 cos 09 (cos 9) fc ■ j a {x) 



[i I i l 2 mlfcl 2 

l4 + Mk \y\ (5 23) 

6(cos9) fc - 9(cos0)p &jo{x) 



mlfcl 



Note that this differential scattering rate can be nonzero 
only on the forward-scattering side again and depends 



4 The situation is different if the \ atom has a mass different from 
that of spin-tr fermions; m x ^ ra a . By considering the two-body 
scattering in which the \ atom with the initial momentum k 
collides with a spin-cr fermion at rest, the energy and momentum 
conservations constrain the final momentum p to be on a sphere 



-fe 



-fc . Therefore, when m Y > 



defined by p ; k i< . unnmir. nni-n (», 

m a , the \ atom can be scattered into an angle range cos 8 > 
^/l — (m CT /m x ) 2 only, while when m x < m CT , it can be scattered 
into any angles. 



on the azimuthal angle tp because of the p ■ j a term. It 
is easy to check that the integration of the differential 
scattering rate in Eq. (|5.20j) or (|5.23j) over the solid angle 
ciO = d cos dip reproduces the total scattering rate in 
Eq. (I5T71) . 



3. Contribution of contact density 

The contribution of the contact density of a spin- 1/2 
Fermi gas to the differential scattering rate of the x atom 
can be extracted in a similar way. As we discussed at the 
end of Sec. IIII El the last term in Eq. (|5.17l) represents 
the contribution from the three-body scattering of the x 
atom with the initial momentum fc with a small pair of 
spin- up and -down fermions at rest. By using the optical 
theorem, the imaginary part of the forward three-body 
scattering amplitude t x (k; fc) can be written as a form of 
the total scattering rate: 



2Im£ x (fc;fc) 

'dpdq^dq^ 



(5.24) 



t x (k;p)A(e k - e p ,k-p) 



(2tt) 9 

+ t t (k;q t )A i (e k -e qf ,k- q t ) 
+ h(k; qi)A t (e k - t qi , k - q ; )| 2 
x (2tt) 4 (5(p + q t + Q| - k)5(e p + e qt + e Ql - e k ). 

Here, p and q a are momenta of the x atom and the spin- 
cr fermion in the final state, respectively. This equality 
can be checked by a direct calculation starting with the 
set of integral equations (j5.10[) (see Appendix [C} . 

The right-hand side of Eq. (|5.24|) is infrared diver- 
gent at \q a \ — > because of t a (k;q a ) — > m/q 2 for both 
a =til [ scc Eq. (15. lOjt ] . These infrared divergences are 
canceled exactly by the second term in Eq. (|5.6j) because 
its imaginary part also contains the same form of infrared 
divergences which can be seen from 



2W t (e fe ,fc) 



dq ( m 



(2tt) 3 \q 2 



x (2n) i 5(p + g t - k)5(e p + e qt - e k ), 



(5.25) 



and the same for j" •<->• J,. Therefore, the imaginary part of 
the regularized on-shell three-body scattering amplitude 



2Im^ og (fc;fc) = 2Imt x (fc;fc) 



dq fm\ 2 (5-26) 



2lm[A t (e k7 k) + A i (e k ,k)} J -0- 



appearing in T x (fc) is free of divergences. By recalling 
that p in Eqs. (I5.24[) and (|5.25[) corresponds to the mo- 
mentum of the x atom in the final state, the contribution 
of the contact density to the differential scattering rate 
of the x atom can be identified as 
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FIG. 9. (Color online) Universal function g(cos9, \k\/n' t ) to determine the contribution of the contact density to the differential 
scattering rate (|5.28[) . g(cos8, \k\/n' t ) is periodic in terms of < (so/tt) ln( | fc | / ac^ ) < 1 and here its dependence on cos# is 
shown at four values of (so/tt) hi(|fc|/fv',,) = 3/16, 7/16, 11/16, 15/16. These four values roughly correspond to the minimum, 
inflection point, maximum, and inflection point of the total scattering rate, respectively [see Eq. (|5.29[1 and Fig. [7]. 



dp 



(2tt) 9 



\t x (k;p)A(e k - e p ,k - p) +t t (k;q t )A l (e k - e qt ,fe - q t ) 



ti(k;qi)A t (e k - e q ^k - q±)\ (2tt) 8(p + q t + q l - k)S(e p + + e qi - e k ) 

2 



C(x) 



- ) \A i (e k ,k)\ z (27rr6(p + q l -k)6(e p + e qi -e k ) 

- (™ ) \A t (e k ,k)\ 2 (2n) i S(p + q t - k)S(e p + e qf - e k ) 



(5.27) 



This differential scattering rate can be evaluated by 
using the on-shell three-body scattering amplitudes 
t x ^,i(k;p) at infinite scattering lengths (a^|fc|) _1 , 
(a,|,|fe|) _1 , (a, | A; | ) 1 = obtained from the numerical 
solutions of tp_B{k;p) in Eqs. (|5.11|) - (l5.13p [note that 
t-f(fc;p) = t±(k;p) when a,f = ajj. The physical mean- 
ing of the subtracted terms was discussed at the end of 
Sec. lllTFl 



Finally, by performing the integration over the magni- 
tude of p, the angle distribution of the scattered x atom 



is given by 



gf(fc) 



(C) 



/TP 

d\ P \\ P r L > 



(k) 



g\ cost 



dp 

mk 2 



(5.28) 



The resulting universal function g(cos 6>, |fc|//i'„) de- 
pends on the polar angle 9 and the momentum to 
Efimov parameter ratio |fc|//t* in a log-periodic way. 
g(cos9, |fc|/«») as a function of cos 9 is shown in Fig. [9] 
at four values of (s /tt) ln(|fe|/<) = 3/16, 7/16, 11/16, 
15/16. Note that g(cos6, |fc|/«£) is mostly negative on 
the forward-scattering side (cos 9 > 0) , while it is positive 
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everywhere on the backward-scattering side (cos 6* < 0). 
The divergences of g(cos9, |fe|/«£) at cos 9 = and 1 sig- 
nal poor convergences of the large-momentum expansion 
around these angles. Higher order corrections would be 
important as well and ideally need to be resummed. How- 
ever, away from these two singularities, we expect our 
large-momentum expansion to be valid in a wide range 
of momentum based on the observation in Sec. Ill CI 

The integration of the differential scattering rate (|5,28|) 
over the solid angle dn — d cos 9 d(p reproduces the total 
scattering rate in Eq. (I5.1(: 



T<£\k)=2ii [ d cos 9 g ( cos 9, ^ 
x 7-i V < 

C(x) 



C(x) 
mk 2 



(5.29) 



mk 2 



which relates the integral of g(cos9, |fc|/«£) to the imag- 
inary part of /(|fc|//t*). Note that the above four val- 
ues of (s /7r)ln(|fej/<) = 3/16, 7/16, 11/16, 15/16 are 
chosen so that they roughly correspond to the mini- 
mum, inflection point, maximum, and inflection point 
of Im/(|fe|/«£), respectively (see Fig. [7]). Therefore, al- 
though differences are difficult to see in the differential 
scattering rate from Fig. [SJ the contribution of the con- 
tact density to the total scattering rate varies by up to 
17% for a different momentum to Efimov parameter ratio 
|fc|/<. By combining the results in Eqs. (j5.20j) . (|5.23[) . 
and (|5.28p . we obtain the differential scattering rate of 
the x atom shot into a spin-1/2 Fermi gas up to 0(k~ 3 ): 



dn 



E 

<r=t,4. 



dr { ™° ] (k) dT^">{k) 



CM/ 



dr { x c \k) 
dn 



dn 



o(k- 3 ). 



dn 



(5.30) 



This result was presented previously in Eq. (|2.6I) . where 
the time dependence of the densities is suppressed by 
assuming a stationary state. 



D. Spinless Bose gas 

The analysis in the case of a spinless Bose gas is similar. 
Again a probe atom is denoted by \ and assumed to 
interact with a boson by a scattering length a x . The 
Lagrangian density describing such a problem is 



+ ^Ax + <bcxHHx + £-B, (5.31) 
2m x J 



where Cb defined in Eq. (|4.1I) describes spinless bosons 
interacting with each other by a scattering length a. For 



5 Since dTy" (fc)/dn has a pole at cos 9 = due 
g{cos9, |fc|/ft'„) ~ — 4/(-7r 2 cos 9), the integral over cosd is 
dcrstood as the Cauchy principal value. 



to 



simplicity, we shall assume that the x atom and bosons 
have the same mass m = m x . In analogy with Eq. (|3.5|) . 
the two-body scattering amplitude between the % atom 
and a boson is given by 



A x (k) = 



47T 



1 



K- - mk - i0+ - 

4 u a v 



(5.32) 



The behavior of the single-particle Green's function of 
the x atom, 



ig x (x;k)= I dy e tky (T[ X {x + %)xHx 



I)]), (5.33) 



in the large-energy-momcntum limit k — > oo can be un- 
derstood by using the operator product expansion as in 
Sees. IIIII and IIVI Since the x atom and two identical 
bosons wit h zero-ran ge interactions suffer from the Efi- 
mov effect jllOl llllj . we only consider local operators 
with Nq — up to scaling dimensions Aq = 4 (see dis- 
cussions m Sec. HVB| . Then in analogy with Eqs. (|4~20)) 
and (|4.21l) . Q x ix\ k) can be written in the usual form 



G x (x;k) = 



1 



ko — £fc — S x (x; k) + ^ 



(5.34) 



where S x (a;; k) is the self-energy of the x atom given by 



S x (x; k) - -A x {k) n(x) + ^A^k) ■ j(x) 



T™s(fc,0;fc,0) 



d_ 

dk' 

C{x) 
m 2 



(5.35) 



Here, n(x) = (ip^ip) and j(x) = (—iip^Vip) are the num- 
ber density and current density of bosons and C(x) = 
m 2 (<jr4>) is the contact density of a spinless Bose gas. 
Note that these parameters only refer to many-body 
properties of the given spinless Bose gas and do not in- 
volve the information related to the x atom. On the 
other hand, T™ g (fc, 0; k, 0) is a finite quantity defined by 



dq 

(2^7 



7^(A:,0;fc,0) = T x (*;,0;A:,0)-^ x (fc) 

(5.36) 

where T x (k,p; k' ,p') is the three-body scattering ampli- 
tude between the x atom and a dimer composed of two 
identical bosons with (k,p) [(k',p')) being their initial (fi- 
nal) energy-momentum (see Fig. [5]) . Corrections to the 



above expression of S x (a;; k) denoted by 



start with 



~ (0}/k A °~ 2 , where O are all possible operators with 
the scaling dimension Aq — 5. 

Then by setting ko — in Eq. (|5.35p . the on-shell 
self-energy of the x atom for an arbitrary state is found 
to be 



Y, x (x;e k ,k) 



4n 



n[x) 



a-y \h\ 



t™Hk;k) 



C(x) 



0{k~ A ). 



2iri k ■ j(x) 
7 m\k\ 2 

(5.37) 



v 2 a x | fe| 
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The quasiparticle energy and scattering rate of the x 
atom in a spinless Bose gas are given by the real and 
imaginary parts of H x {x\ek,k) according to Eqs. (|5.8j) 
and (|5.9[) , respectively. Our next task is to determine 
the regularized on-shell three-body scattering amplitude 
f x e Hk;p) = T^(k,0;p,k-p)\ ko=eh , po=ep in Eq. (JOT), 
which requires solving a three-body problem. 



E. Three-body problem 

We now compute the three-body scattering amplitude 
T x (k, 0; fc, 0). Unlike the previous cases in Sees. IIIIFl and 
HVEl T x (k;p) = T x (k,0;p,k-p) by itself does not solve 
a closed integral equation. To find a closed set of inte- 
gral equations, we need to introduce another three-body 
scattering amplitude T^(k;p), which describes a process 
where the \ atom and a dimer composed of two identical 
bosons with their energy-momentum (fc, 0) are scattered 
into a boson and a dimer composed of the x atom and 
the other boson with their energy-momentum (p,k—p). 
These two scattering amplitudes are solutions to a closed 
set of integral equations depicted in Fig.[8]with the iden- 
tification of T-f , Tj_ with T.jf, . Then by following the same 
procedures as in Eqs. (|3.43|) - (|3.45[) . the integral equa- 
tions solved by the on-shell three-body scattering am- 
plitudes t x ^{k;p) = T x ^(tk,k]€ p ,p) can be written as 



and 



t x {k-p) 
and 



dq 

(2^7 



IC ax (k;p, q) mk;q) 



(5.38a) 



. , lib 

U(k;p) = -J 



dq 



JC a (k;p, q) t x (k;q) 



+ 



dq 



(5.38b) 



ICa x {k;p, q) tij,(k;q), 



where the integral kernel K, a (k;p,q) is defined in 
Eq. (|3.46[) . Note that the factor 2 in front of JC a orig- 
inates from the fact that the two-body scattering ampli- 
tude between two identical bosons is twice as large as that 
between two distinguishable particles [compare Eqs. (I3.5|) 
and imp ]. 

As long as we are interested in t r x s (k;k) ~ m/k 2 up 
to 0(k~ 3 ) [see Eq. (|5.37[) ]. we can set the two scatter- 
ing lengths infinite (a x |fc|) _1 , (a|fc|) _1 = 0, because the 
dependence on them appears only from 0(fc~ 3 ). In this 
case, by defining 



t F (k;p) = 2t x (k;p) - t^{k;p) 



and 



t B (k;p) = t x (k;p) + mk;p), (5.40) 

the two coupled integral equations (|5.38l) can be brought 
into two independent integral equations: 

in f dq 

t F (k;p) = -— - J ^-p-/C 00 (fe;p,q)t F (fc;q) (5.41a) 



771 f do 

t B (k;p) = — + 2 7 ^K 00 (k;p,q)t B (k;q). (5.41b) 
p z J (27r) ci 

Because these two integral equations are equivalent to 
Eqs. (f3T45|) and (|43B]) at (alfeD^ 1 = 0, their solutions 
are already obtained. 

By using the numerical solutions obtained previously 
in Eqs. (|330]) and (|4~30f . t™ s (k;k) within the accuracy 
of 0(k~ 3 ) is found to be 



k) = i [t^(k; k) + t^(k- k)] + 0(k- 3 ) 



3.771 + / 



fc 2 



0(fc" 3 ). 



(5.42) 



Here, f(\k\/n' r ) is the universal log-periodic function 
plotted in Fig. [7] and approximately given by Eq. (|4.31j) 
with being the Efimov parameter associated with 
a three-body system of the x atom with two identi- 
cal bosons. By substituting this numerical solution of 
the three-body problem into the on-shell self-energy in 
Eq. (|5.37l) , the large- momentum expansions of the quasi- 
particle energy and scattering rate of the x atom in a 
spinless Bose gas are obtained from Eqs. (15.8[) and (15.91) : 



E x (x;k) 



1 + 32tt 



n(x) 



2 f (\k\ 
- 3.771 + Re/ (U 



C(x) 



0(k~ 5 ) 



(5.43) 



and 



T x (x;k) 



32tt 



n{x) ^ k-j(x) 
Ifcl 3 



Ifel 



Imf 



k\\ C{x) 

w 



0(k~ 5 ) 



(5.44) 



efc, 



respectively. 

Furthermore, by comparing Eqs. (|5.39|) - (|5.41[) with 



(5.39) Eqs. (l5TTT|) - ([5TT3l . we find t^(k;p) = t t ^(k;p) and 
t x (k;p) in a spinless Bose gas is a half of that in a spin- 



1/2 Fermi gas. Therefore, the contribution of the contact 
density to the differential scattering rate of the x atom 
in a spinless Bose gas is also a half of that in a spin-1/2 
Fermi gas (see also Appendix [Cj. All discussions given 
in Sec. IV CI for a spin-1/2 Fermi gas apply equally to a 
spinless Bose gas with minor modifications. In partic- 
ular, the differential scattering rate of the x atom shot 
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into a spinless Bose gas is given by 



^r x (fc) 



32 cos 6 6 (cos ( 



n{x) 

w 



32 *^2 cos 6 9(cos 6) k - S(cos 9) k 
+ 6(cos6»)p} 



g cosy, — 

Hi 



5{x) 
|fc| 4 

fc|\ C(x) 



\k\ 



0{k- 



2m 



(5.45) 



Compare this result with that for a spin- 1/2 Fermi gas 
in Eq. ([2~6]) . 



VI. WEAK-PROBE LIMIT 
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FIG. 10. Schematic of valid regions of our hard-probe formula 
(|6.1[) and the known weak-probe formula (|6.5[) in the plane of 
the incident momentum k and the scattering length a x ~ 
at ~ a i between the probe atom and an atom constituting 
the target atomic gas. 



A. Spin- 1/2 Fermi gas 

In the previous section, the differential scattering rate 
of a different spin state of atoms shot into an atomic gas 
was derived in the systematic large-momentum expan- 
sion. For a spin-1/2 Fermi gas, it is given by 



dp 



+ 



E 

rfrf(fc) 

dp 



dv^\ k ) , dvr\k) 



dp 



Oik- 



dp 



where the three leading terms were obtained in 
Eqs. (|STTg|) . fST2"2]) . and ([5T2T|) . respectively. On the other 
hand, there is another limit in which a different system- 
atic expansion is possible; that is, the limit of a a — > 
where the probe atom interacts weakly with atoms con- 
stituting the target atomic gas. In this "weak-probe" 
limit, the self-energy of the x atom in a spin-1/2 Fermi 
gas can be expanded perturbatively in terms of a a : 



s x( fc ) = ~ E Cc,niy 
o-=t>4- 

dp S aa > (k - p) 



E c - C 



(2tt) 3 po - e p + i0+ 



+ 0(al). 



(6.2) 



Here, c a = —4ira a /m + 0(a 2 a ) is the coupling strength 
between the x atom and a spin-cr fermion and 

1 



Saa> (k) — 



2n 



dx e ikx 



h a (x)h a >m (6.3) 



is a dynamic structure factor of the spin-1/2 Fermi gas 
with translational symmetries assumed again. As before, 
the imaginary part of the self-energy gives the scattering 
rate of the x atom as 



r x (fc) = -2Im[S x (e fe ,fe)] 



(6.4) 



E 

<T,<T' 



CcrCfj 7 



dp 



Saa'(ek - e P ,fc - p) + 0{a%). 



Because p corresponds to the momentum of the x atom 
in the final state, its differential scattering rate can be 
identified as 



^r x (fc) 

dp 



E 



(2nf 



S aa ,(e k -e p) k-p)+0(al), (6.5) 



which is well known in the context of inelastic n eutron 
scattering in condensed matter physics jll 2L Ill3j . 

This formula (I6.5[) for the differential scattering rate is 
valid in the weak-probe limit a a — > but for an arbitrary 
incident momentum k, while our formula (I6.1|) derived in 
Sec. |V]is valid in the hard-probe limit |fe| — > oo but for 
an arbitrary scattering length a a . Therefore, they cover 
different regions in the plane of a a and \k\, while both the 
results become valid in the double limit of a a — > and 
|fe| — >• oo (see Fig. [T0|). In this section, we show that the 
weak-probe limit of Eq. (|6.1|) is indeed equivalent to the 
hard-probe limit of Eq. (|6.5|) . This serves as a nontrivial 
check of our results presented in the previous section as 
well as a new derivation of the large-energy-momentum 
expansion of the dynamic structure factor S aa > . 



B. Dynamic structure factor 

For simplicity, we shall consider a spin-1/2 Fermi gas at 
infinite scattering length a —¥ oo only. In the weak-probe 
limit a a — > 0, we can make the following expansions in 
each term for the differential scattering rate in Eq. (16.11) : 



O(al), 



t a (k;p) = — + O(a 



(6.6) 



(6.7) 



t x {k;p) 



dq m 



C T + C 1 



(27r) 3 q 2 e p + e q + th-p-q - e* - ?'0+ 



+ 0(al). 



(6.8) 



23 



Accordingly, the contributions of the number density 
(|5.19[) and the current density (|5.22p become 



dp 



(2nf 



S(cj - e K ) + 0(4) 



and 



dp 



(2tt) 2 dK 



(6.9) 



(6.10) 



respectively, where we denoted the energy-momentum 
transfer to the medium by (uj,K) = (efc — e p ,fc — p). 
Similarly, after lengthy but straightforward calculations, 
the contribution of the contact density (|5.27l) becomes 



<flf } (fe) 
dp 



Ci S (SE > (u,K)+ (ct . + g )2 5( AL > (w, K) 



+ 



(2tt) 2 v ' ' (2tt) 2 



(2nf 



(6.11) 



where we defined! 
(u,K) 



e("-¥) 

(2tt) 2 



k: 2 



/ a: 2 



2\2 



2g 2 



t S(uj - e K ) 



.12) 



These three functions correspond to self-energy-, 
Aslamazov-Larkin-, and Maki-Thompson-type contri- 
butions, resp ectively, in a direct diagrammatic calcula- 
tion [saliiH. 

Then by comparing Eqs. (16.9[) . (I6.10p . and (|6.11[) with 
Eq. (|6.5p . we find that each spin component of the dy- 



namic structure factor should have the following system- 
atic expansion in the hard-probe limit K — > oo: 



d 

S aa (u),K) = n a 6(uj - ex) + ja ■ ~Qj£<i(uj - e K ) 
+ [5 (SE) (ui, K) + 5 (AL) (w, K)] mC + 0{R- A ) 
for both a = f, | and 
%(o;,X) = Jf) 

= [S^ (w, K) + S< MT ) (w, If)] mC + 0(if- 4 ). 



(6.15) 



(6.16) 



The sum of all four components S(u, K) = 
2 ff(7 'Siii7'(w,if) coincides with the large-energy- 
momentum expansion of the total dynamic structur e 
factor found previously in Refs. [H, HE HE Ill4 Ill5j |. 
This establishes the equivalence between our hard-probe 
formula in Eq. (|6.ip and the known weak-probe formula 
in Eq. (|6.5[) in the limit a„ — > followed by |fe| — » oo 
where both the results become valid. However, we 
emphasize again that, away from this double limit, the 
two results are independent and cover different regions 
in the plane of a a and \k\ (see Fig. [TU)) . 



{27:Y\K\\lmu-^ 



ir 2 0(e K -uj) 



muj+ \K\^Jmu- \- 

- In 



TOW - 



X 2 



(6.13) 



S( MT >(w,K) = 



(2-K) 2 mu}\K\ 



K 1 



(6.14) 



The second term in <S^ se ^(oj, K) originates from the subtracted 
terms in Eq. 1 15. 271 1. While it does not contribute away from 
the single-particle peak at uj = ejf l82fl . this term is essential 
to compute the differential scattering rate in Eq. 1 16. 171 1 below. 
This is because its nonintegrable singularity at q = \K\/2 cancels 
that of the first term in 5( se '(cj, K) which becomes apparent by 
rewriting Eq. 116.121 1 as 
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x [S(u 



dq- 



2q 2 



ff2 \2 



■ 2e q ) - 5(u) - e K ) 



C. Differential scattering rate 

Finally, we present the angle distribution of the large- 
momentum x atom scattered by a spin-1/2 Fermi gas 
at infinite scattering length a — > oo but in the weak- 
probe limit a a — > 0. The integration of Eq. (|6.5[) with 
the dynamic structure factor obtained in Eqs. (16.151) and 
(I6.16[) over the magnitude of momentum |p| 2 <i|p| yields 



dil 



- {(5(cos6»)fe- 6(cos 6) p) j 

[h\\ (cos 6)C + 0{k- 1 
[h n (cos6)C + 0(k- 1 )] 
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.-1M «t a i 



m 
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Here, /iii(cos0) and h^±(cos6) are universal functions 
plotted in Fig. [TT] and defined by 



/>oo 

A,| (cos 0) = 8 / d\p\ \p\ 2 [5 (SE ) (e k -e p ,k- p) 
Jo 

+ S (AL \e k -e p ,k-p)] 



(6.18) 
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FIG. 11. (Color online) Universal functions h\\{cosff) (solid 
curve) and hfi(cos9) (dashed curve) to determine the con- 
tribution of the contact density to the differential scattering 
rate (|6.17p in the weak-probe limit Of , aj. — > 0. 



and 

/>oo 

%(cos0) = 8 / d\p\\p\ 2 [S {AL) (e k -e p ,k-p) 



(6.19) 



These two functions correspond to the weak-coupling 
limit of the function g introduced in Eq. (|5.28[) : 



lim 

a-f ,04. ->0 ' 
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a a \k\ 

1 + h n (cos9) a t ai k 2 , 



(6.20) 



while g was evaluated in the strong-coupling limit 
a^,aj_ — > 00 in the previous section. Note that the Efi- 
mov effect does not appear in perturbative expansions. 
Again, one can see from Eq. (|6.17p that the number den- 
sity and current density contribute to the forward scatter- 
ing only, while Fig. QT] shows /in (cos 0) > h^±(cos6) > 
at cos 9 < so that the contact density gives the leading 
contribution to the backward scattering. 

The integration of the differential scattering rate (I6.17[) 
over the solid angle dfi — d cos 9 dip yields the total scat- 
tering rate in the weak-probe limit a a — > 0: 



T x {k)=An [l 

<r=t,l 



k\ n„ 



k - jo 



m 



[l.60098C + O(fe- 1 )] 



4 + Q ! 
2m 
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+ [9.25387C + O^k- 1 )] ^ + 0{a 3 a ). 



Note that, since dT x (k)/dfl has a pole at cos 9 = due 
to hu(cos9) ~ — l/(7r 2 cos0), the integral over cos 9 is 
understood as the Cauchy principal value. 



D. Spinless Bose gas 

Similarly, in the case of a spinless Bose gas at infinite 
scattering length a — > 00, the differential scattering rate 
of the x atom in the weak-probe limit a x — > is given by 



tff x (fc) 
dp 



(2n) 



;S(e k - e p ,k-p) + 0{a 
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Here, S(uj, K) is a dynamic structure factor of the spin- 
less Bose gas and has the following systematic expansion 
in the hard-probe limit K — > 00: 

d 

S(u,K) = n5(uj- e K ) +j ■ - e K ) 

UK 

+ [S (SE) [u, K) + 2 S( AL > («, K) + S( MT > (w, K)] mC 
+ 0{K- i ). (6.23) 

Therefore, the contribution of the contact density to the 
dynamic structure factor in a spinless Bose gas is a half 
of that to the total dynamic structure factor in a spin- 
1/2 Fermi gas 0. Accordingly, Eqs. (|6T7]) and (|6T2Tj) 
remain valid by setting a<|- = —¥ a x , replacing n-\ + nj. 
(it + j\) with n (J), and dividing the coefficients of C by 
two. 

VII. CONCLUSIONS 

In this paper, we investigated various properties of an 
energetic atom propagating through strongly interacting 
atomic gases by using systematic large- momentum ex- 
pansions. Our main results are summarized in Sec. Hi] 
and consist of the quasiparticle energy and scattering rate 
of an energetic atom in a spin-1/2 Fermi gas (Sec. IIIII) . 
those in a spinless Bose gas (Sec. IIV[) . and the differen- 
tial scattering rate of a different spin state of atoms shot 
into a spin-1/2 Fermi gas or a spinless Bose gas (Sec, ly)) . 
Furthermore, a connection of our hard-probe formula de- 
rived in Sec. [V] with dynamic structure factors in the 
weak-probe limit was elucidated in Sec. I VII 

Our result on the quasiparticle energy in a spin-1/2 
Fermi gas reasonably agrees with the recent quantum 
Monte Carlo simulation [74j even at a relatively small 
momentum |fc|/fcp > 1.5. This indicates that our large- 
momentum expansions are valid in a wide range of mo- 
mentum. Further analysis of quantum Monte Carlo 
data incorporating our exact large-momentum expan- 
sions may allow us better access to the intriguing pseu- 
dogap physics. Also our result on the rate at which 
the atom is scattered in the medium may be useful to 
better understand multiple-atom loss mechanisms due to 
atom-dimer "dijets" produced by three-body recombina- 
tion events [13, HH . We found that the contact density 
has a negative contribution to the scattering rate in a 
spinless Bose gas. This rather counterintuitively means 
that the energetic boson can escape from the medium 
easier than we naively estimate from a binary collision. 



25 



We also proposed a scattering experiment in which we 
shoot a different spin state of atoms into an atomic gas 
with a large momentum and measure its differential scat- 
tering rate (Fig. [T]). Here a nice interplay between few- 
body physics and many-body physics can be seen: The 
angle distribution of the scattered atom is determined 
by few-body physics and its overall magnitude is set by 
many-body physics. We elucidated that, because the 
number density and current density of the target atomic 
gas contribute to the forward scattering only, its contact 
density gives the leading contribution to the backward 
scattering. Therefore, such an experiment can be used 
to measure the contact density (integrated along a clas- 
sical trajectory of the probe atom) and thus provides 
a new local probe of strongly interacting atomic gases. 
Its intriguing analogy to nuclear physics experiments on 
short-range pair correlations in nuclei (32l . l33j should be 
explored further. Also we found that the differential scat- 
tering rate can depend on the azimuthal angle only by the 
current density of the target atomic gas. Therefore, the 
azimuthal anisotropy in the differential scattering rate 
may be useful to reveal many-body phases accompanied 
by currents. We hope this work serves as a promising 
starting point for future ultracold-atom experiments and 
builds a new bridge between ultracold atoms and nu- 
clear and particle physics from the perspective of "hard 
probes." 



and the two-body scattering amplitude between spin-up 
and -down fcrmions is 



2n 



(A2) 



Here, M = m^+m^ is a total mass and \i = m^m^/(m^ + 
mi) is a reduced mass. Results in the case of spinless 
bosons are obtained by removing spin indices and replac- 
ing the two-body scattering amplitude A(k) with that 
between two identical bosons [see Eq. (|4.4|) ]: 



A(k) = 



8tt 



(A3) 



1. One-body and two-body sectors 

The Wilson coefficients of local operators of type 
in Eqs. ()3.9|) - ()3.12|) are determined by matching the ma- 
trix elements of both sides of 

J dy e ifc » T[^{x + f)V{(z - f )] = ]T W 0i (k)Oi{x) 

(A4) 

with respect to one-body states (ip a (p 1 ) | and \ip„ (p)) . The 
matrix element of the left-hand side is given by 
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Appendix A: Derivation of Wilson coefficients 



iG t {k)iG a (p){2>K)H{p-p') 

+ (1 - S ta ) iG t (k - *=g-)iG a {p)iA<!s + v±f) 

x iG a ( P ')iG t (k + 2=lL) e -i(p-p> . 



(A5) 



On the other hand, the matrix elements of local operators 
of type V4-Vv are 

(Mp'MIMp)) =iG a (p)(27r) 4 S(p~p'), (A6) 



(Mp'MMx)\Mp)) = ^G <y (p)^G a {p l y- ^(p ~ P ' )x , 

(A7) 



(Vv(pO|-#j9iVv(z)K(p)) 

= Pi+Pi iG(T(p - )iG ^ )e -i{ P - P ')a 



(A8) 



Here we show how the Wilson coefficients in 
Eqs. (pn31)HpH2]) and Eqs. are derived. For 

generality, we consider spin- 1/2 fermions with unequal 
masses ^ mj. . The propagator of fermion held t/v in 
the vacuum is given by 



1 



ko - 6ka + ^ 



A; 2 



2m„ 



(Al) 



(Mp')\ - ^MM^MAp)) 

= (pi-p' i )iG <T (p)iG !T {p , )e- i ^-v'>, 



(Mp'^-^d.M^lMp)) 



(A9) 



(AlO) 
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(Vv(p )| - w t (fc) = -zG t (fc) 2 A(fc), (A16) 

d+o' , ( Al1 ) " 

= (Pi ^r^^p^G^e-^-P >*, 



2 5 



W - # t^(fc) = -iG t (fc)^A(fc), (A17) 



•/ / (A12) 



(Vv(pOI#^Vv(z)|Vv(p)) 

Tl^-l-T)' , ( A13 ) n 9 

= iG^j^^e-W)-, ^t^^(fc) = -^G t (fc) 2 — (A19) 



= (Po -Po) iG a (p)iG a {p')e 



k{ kj 



(A20) 



Therefore, the expansion of Eq. (|A5I) up to 0(p 3 ) is re- ,,s _ w ,, n _ _ _ 

produced by choosing their Wilson coefficients as W -id i ^l^)W ~ W -d i (^ s d^i)^ ~ W id t ^l^)W ~ u > 

(A21) 

Wt(k) = iG^(k), (A15) and all W (k) = for a =t 

2. Three-body sector 

The Wilson coefficients of local operators of type <^t<^ in Eqs. (|3.9p - (|3.12[) are determined by matching the matrix 
elements of both sides of Eq. (|A4[) with respect to two-body states (</>(p')l an d \4>{p))- The matrix element of the 
left-hand side is given by 



dye ik y{cl>(p')\T[^(x + ^\(x - *)]|#p)> 

= iG t (k)iD(p)(2ir) 4 6(p - p') ( A22 ) 
+ iG t {k - P^L)tD(p)iT t (k - 2^1, £±?L + 2^iL ; fc + £^L 7 E+iL _ Z^L)iD(p')iG^{k + 2^El) e -<(p-p')». 

Here, -D(fc) = — -A(fc) is the dimer propagator and T-f(k,p; k' ,p') is the three-body scattering amplitude between a 

spin-up fermion and a dimer with (k,p) [(k',p')] being their initial (hnal) energy-momentum (see Fig. [SJ. On the 
other hand, the matrix elements of local operators of type ipj^ipi are 

(^p')\mp)) =iD{p){2*) A 5{p~p'), (A23) 



iD(p)iD{p')e-^-^ i J J2_G t (^f q)G l (q + ^)G i (g - *±f) 

dq ( 2/i 



(A24) 



iD{p)iD{p')e-^'> fjp^ I U r). 
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,2\ 



= iD(p)iD(p')e-^-^ ijj^H G t ( p -¥ - q)G±(q + *#-)G^q ^f) (A2g) 



2 

2\ 



= W^W^e-^-*')* i J j^( Pi - p'^^f q)G i (q + **£)G&q - 
= iD(p)i£>(p')e-<<P-*>(p i -p' t ) JJS_(^ 2 + 0(p% 

(0(^)1 -^A^(a:)|0(p)) 

= iD{p)iD{p')e-^-^ i J j0^( Pi p^qjG^ - q)G^q + ^f)G^{q - ^f ) 
= l D{p) l D(j>')e-^-^ l J -^(pt-p'Jipj-p^G^ q)Gi{q+ 2+f )Gi{q _ E+2L) 

= iD(p)iD(2>>)e-^-P'> i J -^ qaGt (^f - q)G i {q + ^)G±(q - 
-UWDWe^fgp^) (^)W). 

W)\W\^){x)\4>{p)) 

= iD(p)iD(p')e-^-^ l J^_( Po - P ' )G t mf q)G i (q + 2±El)G;(g - ^f) 

= o( P 2 ). 

Because the matrix elements of local operators of type <ft (j> are given by 

Wp'MUixMp)) = iD{p)iD(j>')e- i( P-P>, 

(4>( P ')\ - i^d.^(x)\<t>(p)) =Pi±P±iD(p)iD(p')e- i <*-*>'>, 

Wp')\ - id^^ixMp)) = ^ -p^iDipWp^e-^-^^ 
the expansion of Eq. (|A22[) up to 0(p 2 ) is reproduced by choosing their Wilson coefficients as 

W^(k) = ^G(fc) 2 T t (fc,0;fc,0) - W^(k) J 

" -tldiartiS*) 3 J (2^)3 \ q ) "i^a^W 2m t j (2tt) 3 V Q / 
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W_^Jk) = -iG t (kfJLT t (k, P ; k,p)\ p ^ - W_^(k) ^ / ^ (^) 2 , (A36) 



(A37) 



These results are presented in Eqs. p.l3|) - ()3.22p for a spin-1/2 Fermi gas and in Eqs. (|4.14j) - (|4.18|) for a spinless Bosc 
gas. 



Appendix B: Details of solving integral equations 



Here we discuss details of solving the integral equations in Eq. (|3.45p for a spin-1/2 Fermi gas and in Eq. (|4.26[) for 
a spinless Bose gas. 



1. Spin-1/2 Fermi gas 



For generality, we consider spin-1/2 fermions with unequal masses ^ mj., in which the integral equation (|3.45[) 
is modified to 



J 7^33 £a(k;p, q) t\(k;q), 



where the integral kernel tC a (k;p, q) is given by 

2tt 1 



JC a (k;p, q) 



H y/fjik - q) 2 - 2^ (e fct - e qt ) - i0+ - \ £ k - P - q i + e pt + e qt - e fef - i0^ 
For unequal masses, the quantity we would like to compute (|3. 491) becomes 



dq 



K a {k;k,q) 



An 
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(Bl) 



(B2) 



(B3) 



The second term is a simple integral and has the following analytic expression for its expansion in terms of (a|fc|) : 



dq 



K a (k;k,q) 



An 



(m + 1)V2u+1 + 



" M 
(^+1) 3 



- q 2 



2fi 



csin fe) , (u + l) 3 . 
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I + 



2/i 



(B4) 



Therefore, the nontrivial task is the calculation of the last term in Eq. (|B3j) at (a|fc|) _1 = and its first derivative 
with respect to (<z|fc|) -1 , which requires solving the two-dimensional integral equation (|F31[) numerically. 

For numerical purposes, it is more convenient to shift momentum variables as p (q) — > p (q) 
the center-of-mass frame and introduce a dimensionless function 



k to move to 



k 2 



M+m t 



fe 



This function solves a simpler integral equation 

where the new inhomogeneous term I(p) and integral kernel J a (p, q) are defined by 



(B5) 



(B6) 
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and 



Ja(p, q) = /C a (fe; p + j^k, q + 

47T 



M+m t o ni i » n i „2 _i_ „2 _i_ 2m t „ . „ _ "U ip. _ 



(B8) 



M+m t 



Since J7a(p, 5) now depends on the angle only between p and q, each component of the partial- wave expansion 



(p) = ^ 4 (p)P e (cos 9) (cos 9 = k-p) 



solves an independent one-dimensional integral equation 



(p) = -ZW(p)- / dqjW(p,q)s<f\q), 



where the partial-wave projections of the inhomogeneous term and integral kernel are given by 
2£+l f 1 , „ s _, s 2^+1 f 1 , „ P e {cos9) 



Of _|_ 1 z* 1 0/ -I- 1 /" 
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Here, p = |p|/|fe|, q = |q|/|fe| are dimensionless momenta and the mass ratio is denoted by u = mf/mj,. Note that 
the integrations over cos9 can be done analytically by using Gauss's hypergeometric function: 



1 , a Pi(cos9) 
a cos & — = — — 

_! x + ycos9 x(2£+l) 



1)!! V x) 
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2 ' 2 ' + 2'\x) 



(B13) 



In terms of si^ (p) , the last term in Eq. (|B3[) is written as 



dq 



K a (k;k,q) 



H(k;q) + ^ 



E / dqjP(^- 1 ,q)\s ( t l \q)+l^(q) 



2^ 
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(B14) 



E dqjP(£±,q)lW(q) 



2fi 

fc2' 



In the last line, we used the integral equation for s^\p) in Eq. (|B10[) . Because s^(p) is singular at p - 
and tt^tt , it is better to work on the following function: 



tt -(2u+l) 



*4*> (p) ee S f (p)+lW(p)- dqjP(p,q)lW(q), 

Jo 



(«+l)(2«+l) 



(B15) 



in which the singularity at p 



u 2 -(2tt+l) 
(m+1)(2m+1) 



function solves an integral equation 



is eliminated and the singularity at p = 2u+i 1S wea ker. This new 



00 />oo 
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(B16) 
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and the last term in Eq. (IB3|) is given by its value at p — " + • 



2u+l • 



dq 



K. a (k;k,q) 



1=0 



2fi 



(B17) 



In the case of equal masses u = 1, we numerically solved the integral equation (|B16|) at (ak) 1 ~ for < £ < £ n 
up to £ max = 20. By extrapolating the result to £ max — > oo, we obtain 



J2 S S { e) (|) * 2.335 + 7 -^i + 0{k-^ 



1=0 



alfel 



(B18) 



Since we find that the imaginary part of 0(k°) term is as small as ~ 2 x 10~ 6 and the real part of 0(k^ 1 ) term is as 
small as ~ 4 x 10 ~ 5 , we assume that their actual values are zero. This result combined with Eqs. (|B3|) and (|B4[) is 
presented in Eq. (I3.50|) . 



2. Spinless Bose gas 

In the case of spinless bosons, the quantity we would like to compute is Eq. (|4.29l) : 



)C a (k;k,q) - 
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dq 
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in which t(k;p) solves the integral equation in Eq. (|4.26[) : 



t(k;p) 



P 



dq 



)C a (k;p, q)t(k;q) 
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Since the second term in Eq. (|B19I) is twice as large as that in Eq. (|B3|) . its expansion in terms of (a|fc|) _1 is obtained 
from Eq. (|B4[) with u = 1 by multiplying it by two. In analogy with the case of spin-1/2 fermions, we shift momentum 
variables as p (q) — > p (q) + fc/3 to move to the center-of-mass frame and introduce a dimensionless function 



k 2 

s(p) = — t(k;p 
m 



Each component of its partial- wave expansion 

oo 

s(p) = s ( - e) {p)P e (cos9) {cos9 ee k ■ p) 



(B21) 
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1=0 



solves an independent one-dimensional integral equation 



Ml 



(p)=lW(p)+2 dqjP(p,q)sW(q). (B23) 
Jo 

Here the inhomogeneous term X^\p) and integral kernel J [ a (p,q) are obtained from Eqs. (|B11|) and (IB12[) . respec- 
tively, by setting u = 1. 

Then by defining the better behaving function 

I* oo 

Ss^(p) = s^(p)-l^(p)-2 dqjW(p,q)lW(q), (B24) 

Jo 

with its integral equation 

/•OO />oo />oo 

5sM(p)=4 dq dq'J^(p,q)J^(q,q')l^(q') + 2 dq J<f> (p,q) Ss® (q), (B25) 
Jo Jo Jo 

the last term in Eq. (|B19|) is given by its value at p = 2/3: 



dq 



/C a (fc; k,q) 
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We numerically solved the integral equation (|B25|) at (ak) 1 = for 1 < £ < £ max up to £ max = 20. By extrapolating 
the result to £ max — > oo, we obtain 

oo 

Ss W (§) ~ -2.044 + Oik' 1 ). (B27) 

i=\ 

Since we find that the imaginary part of 0(k°) term is as small as <~ 4 x 10~ 6 , we assume that its actual value is zero. 
On the other ha nd, the £ = channel has to be treated specially due to the Efimov effect. 

As is known jl09l llld| . the integral equation (|B25|) in the I — channel is ill defined without introducing an 
ultraviolet momentum cutoff A: 

/>oo /-oo />A/|fc| 

6sM(p)=4 dq dq'j^(p, q )J^( q , q ')l^( q ') + 2 dqjW(p,q)6 S W(q). (B28) 

Jo Jo Jo 

First we find that 5s^°- ) (2/3) obtained by solving this integral equation is a log-periodic function of \k\/A in the limit 
A/|fe| — > oo, which is approximated by 

S 8 (Q)m~ 30216 -12 20 1 ^O 36 cos(2 So In |fc|/A + 3.9604) - z 1.032 sin(2 So In |fc|/A + 3.9604) 

^ ~ ' ' 1 - 0.08460 sin(2s In \k\/A + 3.9604) ' 1 ' 



with so = 1.00624 being the solution to the transcendental equation 

1. (B30) 



8 sinh(fso) 



y/Zs cosh(f s ) 

Then this artificial cutoff A has to be related to the physical Efimov parameter re* defined in Eq. (|4.27|) . To this end, 
we observe that Eq. (|B28[) without the inhomogeneous term 

fc(t»(p) = 1 f A dq g 2 f'dcose P (B31) 

' " '^2 _*L_ i0 +_i 7-i P 2 + q 2 + Pq cosd-^ -»0+ V ' 

is the Skorniakov-Ter-Martirosia n equ ation to determine the binding energy of three identical bosons by identifying 
/c 2 /3 as the collision energy mE [llfjj. By solving this homogeneous integral equation at infinite scattering length 
s->oo, we find an infinite tower of binding energies given by 



_ — 2irn/so 



A ^ 2 



(n-S-oo), (B32) 



v 5.67865, 

from which we can read off the relationship between A and 

A = 5.67865 x k*. (B33) 

This result combined with Eqs. (1B"4|) , (|B19"|) . (jB27l) . and (jB29l) is presented in Eq. (|4~3"Tj) . which determines the 
universal log-periodic function /(|fc|/re*) introduced in Eq. (|4.30|) . 

Appendix C: Derivation of optical theorem in Eq. (|5.24[) 

Here we derive the optical theorem used in Eq. (|5.24l) by a direct calculation starting with the set of integral equations 
in Eq. (|5.10[) . We only consider the case in which a^ 1 = a^ 1 = a^ 1 = and thus i</,(fe; p) = tf(k; p) = t±(k; p) because 
this case is sufficient for the analysis presented in the text. By using the definitions of tF,B(k;p) in Eqs. (|5.1ip and 
(I5.12j) . the three coupled integral equations in Eq. (|5.10l) can be brought into the two independent integral equations 
in Eq. ([B"T3| : 

t F (k;p) = - - y t F (k- Q )A(e k -e q ,k- g) - - _ - _ . Q+ (Cla) 



and 



t B (k;p) = ^ + lf^ Mk; q) A(e k -e q ,k- q) - - - + - - _ . Q+ . (Clb) 
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Here we used the expression of lC a (k;p, q) in Eq. (13.46[) and A(k) = —D(k) is the two-body scattering amplitude in 
Eq. (|3.5I) at infinite scattering length. Their complex conjugates are 

m f dq 



and 

*<*"> - ? + 2 / wr * K p+t , + J p -,- tl , + m - < C2b > 

By using these complex conjugates, t FtB (k;p) at p — k can be written as 

. , , , m f dp m 

t F {k;k) = - - / ——=t F {k;p)A(e k - e p ,k-p) — 

p 1 J (2Try p z 

TTL I (JfO 

= 2+ \tF(k;p)\ 2 A(e k -e p ,k-p) (C3a) 

'dpdq , , f 



Jj7^ t F{k;p)A(e k - e p ,k-p) t* F (k;q)A*(e k -e q ,k- q) 



and 



, , . m f dp m 
t B (k;k) = — + 2 / ——t B (k;p)A(e k - e p ,k-p)^ 
p z J (2ny p z 

= ^ + 2 J-^\tB(k;p)\ 2 A(e k -e p ,k-p) (C3b) 

- 4 / i^^tB(k;p)A{e k ^e p ,k-p)t B (k;q)A*(e k -e q ,k-q) — — — . 

Then by using the identity 

2 ImA(e fe - e p , k - p) = J ^0^\A(e k -e p ,k- p)\ 2 (2ir) 4 6(p + g t + 9i - fe)5(e p + e 9f + e„ - e fe ), (C4) 

the imaginary parts of t FjB (k] k) become 

2Imt F (fc;fc) = i [ dP f q \ d g qi \t F (k;p)A(e k -e p ,k-p)- t F (k; q t )A(e k - e q k - q t )\ 2 

I J (27TJ 3 (C5a) 

x (2tt) 4 (5(p + 9t + 94 - fe)^(e P + e qt + e q+ - efe) 

and 

2Imi B (fc;fc) = \ [ ^f^^ 1 \t B (k; p)A(e k - e p ,k-p) + t B (k; q t )A(e k - e q ,k- q t ) 

3 ^ ( 27r ) (C5b) 

+ t B (k; qi)A(e k - e^, k - | 2 (2tt) 4 (5(p + q t + - k)S(e p + e qr + e qi - e k ). 

Finally, by using the definitions of t FB (k;p) in Eqs. (|5.11l) and (I5.12[) again, the imaginary part of the forward 
three-body scattering amplitude t x (fe; fc) is found to be 



2lmt x (k;k) = - lm[t F (k;k) +t B (k;k)} 



4 
3 

= / P (2^ \ t x(k-,P)A(e k -e p ,k-p)+ fy(fc; q^)A{e k - e qf , k - q t ) (C6) 
+ H{k;q i )A{e k - e qi ,k - q;)| 2 (2ir) 4 6(p + q f + q l - k)S(e p + € qi; + e qi - e k ), 



which is equivalent to Eq. (|5.24|) when a* 1 = a , 1 = a 1 = 
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On the other hand, in the case of spinless bosons, the definitions of tp.B(k;p) in Eqs. (|5.39[) and (15.40[) lead to 

2 
3 

dpdq t dq l 1 , )A( _,u_ v ) + t ,(u. „+)A(e k - e„„. fc - Of) ( C7 ) 



2Imi x (fc;fc) = -lm[t F (k;k) +t B (k;k)} 



= J ^^n)^ \ l 2< x( fc ;P)^( e fe - £ p> k ~ P) + 9t)^( e fc ~£q t ,k- g t ) 
+ iv( fe !9-l-)^( e fe _ e gi' k - 1l)\ 2 (2tt) 4 (5(p + q t + - fe)<5(e p + e qt + e g+ - e fc ). 

Since t x (k;p) in Eq. (|C7|) is a half of that in Eq. (|C6|) . the integrand in Eq. (|C7[) is also a half of that in Eq. (|CGD . 
This establishes that the contribution of the contact density to the differential scattering rate of the x atom in a 
spinless Bose gas is a half of that in a spin-1/2 Fermi gas. 
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